Obscuration in extremely luminous quasars 
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ABSTRACT 

The spectral energy distributions and infrared spectra of a remarkable sample of obscured 



AGN s selected in the mid-infrared are modeled with recent clumpy torus models (jHonig et al 



20061 ) to investigate the nature of the sources, the properties of the obscuring matter, and depen- 
dencies on luminosity. The sample contains 21 obscured AGNs at z=1.3-3 that were discovered 
in the largest Spitzer surveys (SWIRE, NDWFS, & FLS) by means of their extremely red in- 
frared to optical colors. All selected sources show the 9.7/xm silicate feature in absorption and 
have extreme mid-infrared luminosities (i/L,y(6/im) = L(6/im) ~ IC*^ ergss"^). 

The infrared SEDs and spectra of 12 sources are well reproduced with a simple torus model, 
while the remaining 9 sources require foreground extinction from a cold dust component to 
reproduce both the depth of the silicate feature and the near-infrared emission from hot dust. 
The best-fit torus models show a broad range of inclinations, with no preference for the edge-on 
torus expected in obscured AGNs. 

Based on the unobscured QSO mid- infrared MIR luminosity function ( Brown et al.ll200^ . and 
on a color-selected sample of obscured and unobscured infrared sources, we estimate the surface 
densities of obscured and unobscured QSOs at L(6/im)>10^^Lo and 2=1.3-3.0 to be about 17- 
22 deg~^, and 11.7 deg~^, respectively. Overall we find that ^35-41% of luminous QSOs are 
unobscured, 37-40% are obscured by the torus, and 23-25% are obscured by a cold absorber 
detached from the torus. These fractions constrain the torus half opening angle (measured from 
the torus axis) to be ^67° . This value is significantly larger than found for FIR selected samples 
of AGN at lower luminosity (^^46°), supporting the receding torus scenario. 

A far-infrared component is observed in 8 objects. The estimated far-infrared luminosities 
associated with this component all exceed 3.3xlO^^Lo, implying SFRs of 600-3000 M© yr"-'^. 
For the whole sample, the average contribution from a starburst component to the bolometric 
luminosity, as estimated from the PAH 7.7fim luminosity in the composite IR spectra, is <20% 
of the total bolometric luminosity. 

Five sources are detected in the X-rays, and large column densities (Nh>10^'^ cm^^) are 
estimated for all of them. 

Subject headings: galaxies; active — - galaxies: high-redshift — quasars: general — infrared: galaxies 
1. Introduction 
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1.1. Obscured AGNs at high luminosities 

Multiple X-ray studies show that the frac- 
tion of active galactic nuclei (AGNs) whose emis- 
sion is heavily absorbed (Nh> lO^^cm"^) de- 



increasing luminosity (from >80 
ergss-i, to 38% at lO^^ergss-^ 
Akvlas et al.li2006,l . However, because of the difh- 
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culty in detecting and identifying absorbed AGNs, 
it is still unclear whether their paucity at high 
luminosities is an observ ationa l selection effect 

20051: 



or real ( La Franca et al 



Treister &: Urrv 



2006t lAkvlas et all l200d [ibzzi et alj |2006ll In 

order to overcome these uncertainties, searches 
for absorbed QSOs have focused on observations 
at wavelengths less affected by absorption, i.e., 
infrared (IR) and radio, e.g., FIRST, 2MASS and 



various Spitzer surveys (e.g. iPolletta et al. 
Martmez-Sansigre et al. 2006 : Lacv et al.l 



2006; 



2007 



Wilkes ct al. 2002; Ur rutia et al.ll2005l ). 

These searches have unveiled a large population 
of QSOs obscured at optical wavelengths. As- 
suming that all optically obscured QSOs are also 
absorbed in X-rays, the frac tion of absorbed QSOs 
woul d be >50% of aU QSOs (iMarti'ncz-S ansig re et al 
20051 ). This fraction is still significantly lower 



than the fraction m easured for AGNs at lower 
luminosities, 80% (jOsterbrock &: Martell Il993l: 
lAkylas et al. 200^). 

Current AGN evolutionary models and observa- 
tions indicate a link between the absorption prop- 
erties in AGNs and their lumin osities (see e.g. Fig- 
ure 4 in iHopkins et al.l 120051 ). More specifically, 
AGNs in a growing phase are moderately lumi- 
nous and heavily absorbed, and as the AGN lu- 
minosity increases the intense radiation field de- 
stroys the surrou nding matter and the A GN shines 
unabsorbed fe.g. iDi Matteo et al.l 120051 ). A simi- 
lar tren d is predicted by the receding torus mod- 
els (e.g. Lawrencell991 : SimpsonlboOS ). According 
to these models, the opening angle of the torus 
(measured from the torus axis to the equatorial 
plane) is larger in more luminous objects. Thus, 
obscuration should be less common in more lumi- 
nous AGNs. 

1.2. The obscuring matter in AGNs 

Obscuration in AGNs is caused by a mixture 
of neutral and ionized g well as dust. Ab- 

sorption by gas is usually observed in the X- 
rays through the suppression of soft X-ray emis- 
sion with an energy cut-off that increases with 
the gas column density, Nh- The absorbing gas 
is mainly located in the circumnuc lear region as 



suggested b y measured variability (jRisaliti et al 



2005L |2003), but gas in our Galaxy and in the 
host galaxy can also contribute to the overall ab- 
sorption. Obscuration by dust is usually observed 



at optical and IR wavelengths where the contin- 
uum and broad emission lines are reddened and, 
in some cases, completely suppressed. In the near- 
and mid-IR (NIR and MIR), dust obscuration in- 
terplays with thermal re-emission from the puta- 
tive parsec-scaled dust torus. Depending on the 
actual geometry, optical depth and chemistry, dust 
absorption features at about 10 and 18/im due to 
silicates (Si) can be present in the IR. 

According to AGN unification mod els (e.g. 

'Antonucci"l993') and to torus mode ls (e.g., Krolik &: Begelman 
1988; Pier fc Krolik 1992: Efstathiou fc Rowan-Robinsonl 
19951 : iGranato et all Il997l : 



Nenkova et al 



DuUem ond fc van Bemmel 2005t Fritz et al 



2002; 



20061: 



Honig et al. 2006f ) optical depths are expected 



to be at least roughly correlated from the X-ray 
through the MIR, although this correlation might 
be very weak for clumpy tori. In obscured AGNs, 
the line of sight (LOS) intercepts the putative 
torus and thus absorbs radiation from the nucleus 
(broad emission lines. X-rays and thermal radia- 
tion from hot dust), while in unobscured AGNs the 
radiation from the nucleus is directly visible and 
thus unabsorbed. An edge-on view of the torus 
is characterized by red optical through IR colors 
and Si features in absorption (at 10 and 18/im), 
while a face-on torus has bluer colors and Si fea- 
tures absent or in emission. Sources with narrow 
emission lines in their optical spectra, red colors 
and strong Si absorption features are usually as- 
sociated with a highly inclined torus and thus are 
expected to be also absorbed in the X-rays. How- 
ever, there is emerging evidence for a surprising 
mismatch between the absorption measured in the 
X-rays and that measured in the optical through 
IR for a significant fraction of sources, especially 
at high luminosities. 

AGNs selected in the mid-IR because of ex- 
tremely red optical through IR colors sometimes 
show broad optical emission lines, even hav- 
ing absorpti on features in the ir IR spectra in 
some cases (jBrand et al.l 120071) . The compar- 



ison between X-ray absorption and optical ob- 
scuration in various X-ray selected AGN samples 
shows that about 20-30% of obscured AGNs are 
not absorbed in the X-rays and vi c eversa (e.g. 



IPcrola et all 12004 



t2007; Ghozzi et al 



Tozzi et all l2006t iTaier et al 



20071) . The MIR spectra of a 



sample of type 2 QSOs with heavy X-ray absorp- 
tion did not reveal the Si feature in absorption 
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at 10/im ( Sturm et al. 20061 ) as expected. This 
mismatch is even more pronounced at high lumi- 
nosities. Indeed, absorption signatures in the IR, 
e.g. the Si absorption feature at 10/^m, are more 
prominent at the high luminosities of ultra lumi- 
nous IR galaxies (ULIRGS) (San ders fc Mirabel 
Il996l : iSao et al.ll2007h . while X -ray absorption is 
progressively les s common in AGNs with increas- 



ing l uminosity (jUeda et al.l l2003t iHasinger et al 



20051 ). These results suggest that a non negli- 
gible fraction of obscured AGNs might not be 
obscured by a to rus, but by dust in either the nar- 
row line region (jSturm et al.ll2005( ) or in the host 
galaxy and that large X-ray column densities are 
not always associated with geometric ally and opti- 



cally t hick dust distributions (see alsolRigbv et al 



I2OO6I: iMartfnez-Sansigre et al] I2OO6I ). Alterna- 
tively, this mismatch can be explained by a low 
dust-to-gas ratio (Ay/Nn) or by a different path 



for the IR and the X-ray LOSs (jMaiolino et al 
200ll:IShi et alll2006l) . 



Large optical depths for Si absorption features 
(>1.7) usually imply a compact source deeply em- 
bedded in a smooth distribution of material which 
is both geometrically and optically thick, rather 
than absorption by a foreg round screen of dust 
close to the heating source (jLevenson et al.ll2007 : 
Imanishi et al. [ioOTh . Also a detached foreground 



cold absorber which is far enough from the heat- 
ing source to be cold can produce large Si optical 
depths. Deep {rsi >1.7) Si absorption features 
are usually considered as indicators of a buried 
compact source such as an AGN. However, it is not 
clear whether the optical depth is a consequence of 
a random alignment of Si clouds or a specific dust 
distribution and orientation. The observed MIR 
spectra of obscured AGNs do not favor the latter 
scenario because the Si absorption features are not 
always present in the spectra of optically obscured 
and X-ray absorbed AGNs, as would be expected 
if the Si absorbers were associated with the same 
material which suppresses the broad emission lines 
and the soft X-ray emission. 

In order to investigate the properties of the ob- 
scuring matter in extremely luminous AGNs and 
asses how often the absorption signatures at op- 
tical, X-ray and infrared wavelengths do not cor- 
relate, a comprehensive analysis of the absorption 
properties at all wavelengths of obscured AGNs 
at high luminosity is necessary, as well as detailed 



studies of their absorbing matter. Observations 
at various wavelengths from X-ray through opti- 
cal and IR constrain the absorption along multiple 
LOSs and thus the geometry and distribution of 
the absorbing matter. 

In this work, we investigate the observed prop- 
erties of the obscuring matter in a sample of 
obscured and extremely luminous AGNs, and 
we model their SED s with clumpy torus mod- 
els ( Honig et al. 20061 ) to explore the dust geome- 
tries associated with large obscuration. This study 
is based on a sample of extremely luminous and 
obscured AGNs for which multi-wavelength SEDs 
and IR spectra are available. The selected sample 
includes sources from the three widest Spitzer ex- 
tragalactic surveys which were obser ved with the 



Spitz er InfraRed Spectrograph (IRS; iHouck et al 
20041 ) ■ In § [21 we describe the sample selection, 



and in § |3] the observations and the data used in 
this work. The IRS spectra and detected spec- 
tral features are presented in § |4l The modeling of 



the S EDs using clumpy torus models (jHonig et al 



20061 ) is presented in § [S] The composite spectra 



of two sub-samples defined on the basis of the 
model results are analyzed in § [B] The proper- 
ties (IR colors, surface density, luminosity, optical 
depths, and redshift distribution) of the selected 
sample are compared with those of other samples 
of AGNs and ULIRGs from the hterature in § [71 
The X-ray properties and a comparison between 
the absorption seen in the X-ray and the optical 
depth in the IR are described in § [51 Our results 
and their implications are discussed in § [9l and 
summarized in § 1101 

Throughout the paper, the term "absorbed" 
refers to X-ray sources with effective column den- 
sities NH>10^^cm~^, and "obscured" to sources 
with red optical-IR colors, e.g. F(3.6/im)/F(r')>20 
or F(24/^m)/F(r')>100, implying extreme dust 
extinction. The IR SEDs are defined as "red" 
if they are as red or redder than a power-law 
model, Fi, oc with slope a=2. The terms 

type 1 and type 2 refer to AGNs with broad 
and narrow optical emission lines, respectively, 
in their optical spectra. We adopt a flat cosmol- 
ogy with Hn = 71kms~^Mp c~\ flM=0.27 and 



fiA=0.73 (jSpergel et al.ll200: 
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2. Sample description 



There are three large area Spit zer surveys 
with Infrared Array Camera (IRAC; iFazio et alj 
2004). Multiban d Imaging Photometer (MIPS; 
Rieke et al. and optical photometric cover- 

age, the 50deg^ Spitzer Wi de- Area Infrared Ex- 
tragalactic Survey (SWIReS; iLonsdale et aDl2nn,l 
20041) : the 9 deg^ NO AO Deep Wide-Field Sur- 
vey fNDWFS) of the Bootes field (' jannuzi fc Dev 
IQQalMurrav et al. 2005): and the 3. 7 deg^ Extra- 



galac tic First Look Survey (E-FL^; iFadda et"al 
20061 ). 

The initial objective for Spitzer IRS spec- 
troscopy of sources chosen from these surveys was 
to understand populations of new sources that 
would not have been known prior to Spitzer. Con- 
sequently, spectroscopy has emphasized sources 
which are optically very faint, typically mag(i?)>24, 
and IR bright, i.e. F(24/im)>l mJy. Various se- 
lection criteria based on IRAC colors also entered 
choices of spectroscopic targets, but the unify- 
ing theme of existing spectroscopic samples is 
their high ratio of mid-infrared to optical flux, 
IR/opt = i^Fi,(24^m)/j/F^(i?). Most sources cho- 
sen for spectroscopy have IR /opt > 10. For 58 
sources 



observed in B ootes 
20063). 



(|Houck et al.l l2005t 



Weedman et all l2006cD. 70 sources observed in 
the FLS (|Yan et al.ll2007l : I Weedman et al.ll2006bl ). 
and 20 sources in the SWI RE Lockman Hole 
field (jWeedman et al.l l2006al) . the majority are 



at high redshift (z ^2) and many have infrared 
spectra dominated by the strong silicate {Si) ab- 
sorption feature at rest frame 9.7 /im. The pres- 
ence of this absorption combined with the red 
colors and high luminosities of these sources led 
these previous observers to the conclusion that 
IR bright sources with large IR/optical flux ra- 
tios are dominated by obscured AGN (or QSOs 
because of their very high luminosities). A minor- 
ity of sources in these studies show strong PAH 
emission features characteristic of starbursts, but 
these sources generally do not show strong Si 
absorption, so different populations of AGN and 
starbursts are observed. 

If the interpretation of the obscured sources as 
AGN is correct, this implies a significant popu- 



^http: / /swire. ipac.caltech.edu/swire/swire. html 
^http: / /ssc. spitzer. caltech.edu/fis/ 



lation of luminous, obscured QSOs. It is essen- 
tial, therefore, to verify the AGN classification and 
to understand the obscured sources in context of 
models for AGN. 

Our goal in this work is to determine quantita- 
tive AGN models for the most luminous of these 
obscured sources, with two primary objectives. 
The first is to investigate the properties of the 
AGN emission and the nature of the obscuration 
in AGNs at high luminosities. The second is to 
use the SEDs to define the color criteria which de- 
scribe such obscured AGNs in the IRAC and MIPS 
color space in order to compare surface densities 
between obscured and unobscurcd QSOs at high 
luminosities and high redshifts. To achieve these 
goals, we have selected for detailed modeling the 
most luminous of the obscured sources with AGN- 
dominated MIR SEDs from the archival Spitzer 
spectroscopy in the survey fields mentioned. 

Our selection criteria are the following: high 
redshift, i.e. z>l; large MIR luminosities, i.e. 
6/Ltm luminosities, i/Ly(6/xm) — L(6^m) >10^^Lo; 
the presence of Si in absorption in the IRS spec- 
tra; and red MIR SEDs c onsistent with being 



AGN-dominated (see e.g . IPolletta et all 120061: 



Alonso-Herrero et al.ll2006l ). There are 21 sources 



that s atisfy all these crite r ia; 13 are from the liter- 



ature ( Houck et al.l 120051 : IWeedman et al. l2006d : 



Yan et al. 2007t ). and 8 are from our own IRS ob- 



serving programs. As a consequence of these selec- 
tion criteria, all sources have 24/^m fluxes greater 
than ImJy and have relatively large IR/optical 
flux ratios, i.e. F(24^m)/F(r')>400. To our 
knowledge, this sample contains all sources with 
available IRS spectra in the literature or the 
Spitzer public archive that satisfy the selection cri- 
teria described. However, this sample is not meant 
to be complete or unbiased since it is mainly based 
on the availability of an IRS spectra, therefore it 
is not essential to be exhaustive. 

The selected sample includes the most luminous 
obscured QSOs currently known. The large MIR 
luminosities of the selected targets enable us to 
measure the emission from the hottest dust in the 
torus with negligible contribution from the host 
galaxy, providing a laboratory to test the pre- 
dictions of torus models in obscured AGNs. Up 
to now, these kinds of SEDs could be sampled 
only using high spatial re solution data in ne arby 
AGNs (e.g. in NGC 1068;[H6ni£etaD[200i). By 
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combining samples from the widest Spitzer sur- 
veys, we were able to find a significant number of 
these rare objects. 

The list of selected sources, coordinates, opti- 
cal and IR fluxes are reported in Table [TJ For 
simplicity, we use throughout the paper simplified 
names for the selected sources. Official lAU names 
are reported in Table [TJ For the SWIRE sources, 
the first two letters in their names designate the 
field where they were discovered (LH for Lockman 
Hole, Nf for ELAIS-Nl, and N2 for ELAIS-N2). 
The first letter of the NDWFS source names is B 
for Bootes. The names of the FLS sources corre- 



IRS data were in part taken fro m the archive 



and i n part from our own p rojects (jHouck et al 



spond to those in lYan et al.l (|2007n . 



3. Observations and Data Analysis 

Multi-band photometric data in the optical 
from the ground and in the IR from Spitzer are 
available for all sources. The available data and 
references are listed in Table [T] All sources are 
detected with MIPS at 24 /im with fluxes rang- 
ing from f .0 to 10.6 mJy, and median flux of 2.6 
mJy. All sources have been observed with IRAC 
in 4 bands, from 3.6 to 8.0 /zm with varying lim- 
its. The FLS sources have few detections with 
IRAC, probably because of their less sensitive ob- 
servations compared to the NDWFS and SWIRE 
surveys. Eight sources have been observed in 5 
optical bands, from U- to z-band, 3 in 4 bands, 5 
in 3 bands, and 5 in 1 band. Optical magnitudes 
for Bootes sources are in NDWFS Data Release 3 
(DR3). Optical an d IR fluxes for th e FLS sources 



were ob tained from lYan et al.l (|2007n : ISajina et al 



(|2007ah . and those for the SWIRE sources were 
taken from the latest internal catalogs that will 
be released in the Dat a Release 6 (for d etails on 
the data reduction see Surace et al. I l2005h . 



The IRAC fluxes of the sources in the NDWFS 
survey were measured from the post-BCD images 
in the Spitzer archive. Aperture fluxes were mea- 
sured within a 2'.'9 radius aperture at their position 
and background subtracted. Aperture corrections 
as derived from the SWIRE survey were applied. 
The aperture corrections are 1.15, 1.15, 1.25, and 
1.43, for the four IRAC bands, respectively. 

All sources have been observed in the optical 
r' or R band and 16/21 have been detected. The 
F(24/im) /F(R) flux ratio ranges from 575 to 25000 
with a median value of 3500. 



2005 : Weedman et al. 2006a[ ). The spectroscopic 



observations were made with the IRS Short Low 
module in order 1 (SLl) and with the Long Low 
module in orders 1 and 2 (LLl and LL2), described 
in iHouck et al.l (j2004H ^ providing low resolution 
spectral coverage from ^^8 fim to '^35 /im. The 
main parameters of the IRS observations are listed 
in Tabled! 

Since the objects have similar properties in 
terms of SED and fluxes, we applied the same 
reduction procedure to all spectra for unifor- 
mi ty. The reduction method is described in detail 



in IWeedman et al.l (j2006a[ ) and briefly summa- 
rized here. Six SWIRE sources (see details in 
Table [2]) were processed with version 11.0 of the 
SSC pipeline; the remaining sources were pro- 
cessed with version 13.0. Extraction of source 
spect ra was done w ith the SMART analysis pack- 
age (jHigdon et al. [2004) . Since the selected 
sources are faint and the spectra are dominated 
by background signal, we restricted the number 
of pixels used to define the source spectrum to a 
width of only 4 pixels (scaling with wavelength) 
in order to improve the spectral signal-to-noise 
(S/N). The flux correction necessary to change 
the fluxes obtained with the narrow extraction to 
the fluxes that would be measured with the stan- 
dard extraction provided for the basic calibrated 
data was measured by extracting an unresolved 
source of high S/N with both techniques, and 
this correction (typically ^ 10% but varying with 
wavelength) was applied to all sources. 

A search in the literature provided additional 
data for 6 sources in the sa mple, an ISO 15Atm flux 
measurement for N2_08 ( Gonzalez- Solares et al 



20051 ). and MAMBO 1.2mm flux upp er hmits for 



afl 5 E-FLS sources (|Lutz et al.ll2005h . 

The SEDs and IRS spectra of all sources in the 
sample are shown in Figure [T] and [H The dis- 
played spectra have been boxcar-smoothed to nine 
times the approximate resolution of the different 
IRS modules (0.6 /im for SLl, 0.9 iim for LL2, and 
1.2 /im for LLl). The MIPS 24 /im fluxes typically 



•^The IRS was a collaborative venture between Cornell Uni- 
versity and Ball Aerospace Corporation funded by NASA 
through the Jet Propulsion Laboratory and the Ames Re- 
search Center 



5 




Fig. 1. — Rest-frame spectral energy distributions in uLi, vs A (black symbols and curve) and best-fit model 
(solid grey curve) obtained from a torus model (dashed grey curve) and an elliptical template (dotted grey 
curve) to fit the residuals in the ne ar-IR. Open diamonds represent ISO IS^tm ( in N2_08) and MAMBO 



1.2 mm (in all E-FLS sources) data ( Gonzalez-Solares et al. 20051 : Lutz et al. 2005 ). Downward arrows rep- 



resent 5(7 upper limits for the optical and infrared data points, and 3ct upper limits for millimeter data. 
Source names and redshifts are annotated. 



agree within 10 % of the 24 fim IRS flux. 
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Fig. 1. — Continued 
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Fig. 2. — Observed IRS spectra (grey solid curve) and broad-band IR photometric data (black full circles) 
of all selected sources. The expected location of some spectral features are annotated at the corresponding 
observed wavelength (P refers to PAH emission features). The IRS spectra of the Seyfert 1 Mrk 231 (dotted 
curve) and of the heavily obscured AGN and ULIRGIRAS 000183—7111 (dashed curve) are overplotted for 
comparison. The name and redshift of each source are annotated and rest-frame wavelengths are reported 
on the upper horizontal axis. 



4. Spectral properties of the IRS spectra 

The IRS spectra of all 21 sources are shown 
in Figure [2] The spectra of an abs orbed Seyfert 
1 (Mrk 231; IWeedman eral]|2005f) . and a heav- 
ily obscured ULIRG (IRAS F00183-7111; 100183 



hereinafter, iTran et al.ll200ll; ISpoon et all 120041 ) 

are also shown for comparison. Both refer- 
ence spectra exhibit the Si absorption feature 
at 9.7 ^m but with different optical d epths, 



T5,(Mrk 231) = 0.65 (ISpoon et al.l l2002t ) and 



Ts,(I00183)>5 (jSpoon et al.ll2004f ). The spectrum 
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of 100183 also shows absorption features at 4.26 
and 4.67/im due to CO2 and CO gas, and at 6.85 
and 7.25 fiia due to hydr ogenate d amorphous car- 
bons (HAC) (|Spoon et aLll2004h . 

Spectroscopic redshifts are derived from the 
IRS spectra based on the location of the Si 
feature in comparison with the spectra of Mrk 
231 and 100183. Only in 2 cases (LH_A4, and 
B3) are optical spectroscopic redshifts avail- 
able (|Polletta et al.ll2006t lOesai et al.ll2006l ). For 
most of the sources, IR spectro scopic redshifts 
were previously reported (.Weedrnan et al. I l2n06al : 



Yan et al.ll2007t iHouck et al.ll2005l ). although we 
have revised some of them. The estimated red- 
shifts are reported in Tabled Based on the com- 
parison with previous z determinations and with 
the spectroscopic redshifts, we estimate that un- 
certainties associated with these redshifts are as 
much as ±0.2 for sources with the poorest S/N. 

The Si feature in absorption is clearly observed 
in 18 objects (LH.02, LH.A4, LH.A5, LH_A6, 
LH.A8, LH_A11, Nl_09, N2.08, N2.09, Bl, B2, 
B3, B4, B5, MIPS42, MIPS78, MIPS22204, 
MIPS22303), and only marginally observed in 3 
objects (LH_01, N2_G6, MIPS15840). The ap- 
parent optical depth of the Si feature, Tsi, is 



measured as ln{F^^j^^ 



j rpobs \ 



where Fc 



obs 
9.7 



IS 

the observed flux at 9.7/xm, at the minimum of 
the Si feature, and Fg'^Yfim is the intrinsic flux at 
9.7 Ilia that would arise from an extrapolated, un- 
absorbed continuum. Since in most of the cases, 
the SEDs and spectra do not sample the wave- 
length region beyond the Si feature, we cannot 
use the standard meth od to derive the i ntrinsic 
continuum (see e.g. ISpoon et al. 2004( l. The 
intrinsic flux at 9.7 /im is thus estimated using 
two different methods, 1) by extrapolating a 4- 
8/im rest-frame power-law fit to 9.7 /xm, and 2) 
by normalizing an unobscured AGN template to 
the observed 7.4/im flux, after redshifting it to 
the redshift of each source, and reddening i t using 



the Galactic Center extinction curve (iChiar et al 



20061 ) until the observed Si feature is well repro- 
duced. The Si optical depth, rsi, is measured 
only in sources at z < 2.7, 18 sources. For the 
remaining we estimate a lower limit since we have 
only an upper limit to the flux at the bottom of 
the Si feature. The measured values, reported in 
Tabled vary from >0.2 to 3.4. 

The expected location of other spectral features 
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common in starburst galaxies and AGNs are also 
annotated in Figure [2j 1) absorption features due 
to molecular gas, C02A4.26^m, and A4.67^m, to 
water ice at A6.15 and 7.67/zm, and to HACs at 
A6.85,and 7.25/im; 2) emission features associated 
with polycyclic aromatic hydrocarbons (PAHs) 
at A=3.3, 6.22, 7.7, 8.6, 11.3 ^m, and with low 
and high ionization emission lines, [ArII]A6.99/im, 
H25(3)A9.66/im, and [SIV]A10.54^m. Absorption 
features like CO and HACs would imply the pres- 
ence of warm gas along the LOS and features like 
CO2, water ice and Si would im ply the presence of 



2004 



shield ed cold molecular clouds (jSpoon et al. 
20071 ). Although there are sometimes features in 



the individual spectra that are consistent with 
these spectral features, higher S/N spectra would 
be necessary to confirm their detection. No spec- 
trum shows any secure feature except the Si ab- 
sorption. 

5. Modeling the Spectral energy distribu- 
tions with torus models 

The optical-IR SEDs and re-binned IRS spec- 
tra of all 21 sources are shown in Figure [1] The IR 
SEDs, from 3.6 to 24/^m, and IRS spec tra are mod- 
eled w ith a grid of torus models from lHonig et al 



(l2006h and the residuals in the near-IR (NIR) are 
fit with a galaxy template to represent the host 
galaxy. The model parameters are the cloud den- 
sity distribution (from compact to extended), the 
vertical radial distribution (with various degrees 
of flaring or non-flaring), the number of clouds 
along the LOS ('^optical depth), and the torus 
inclination (from face-on to edge-on). The best 
model is the one that gives the best fit, based on a 
X^-test, to the five broad-band photometric data 
points from 3.6 to 24/zm and to a maximum of six 
additional data points derived by interpolating the 
IRS spectrum at A= 3, 5, 7, 8.5, 9.7, and 12 /xm in 
the rest-frame. 

The best- fit models are shown in Figure [1] and 
the best-fit parameters are listed in Tabled In or- 
der to fit the stellar component, we ad opt a 3 Gyr 



old e lliptical template (from GRASIL; ISilva et al 



19981 ) normalized at the observed flux of the resid- 
uals (model-subtracted SED) at 3.6 /im (in the 
Z-band for LH_A5). The choice of an elliptical 
template is justified by the evidence of ellipti- 
cals as hosts in the vast majority of quasars (e.g. 



Dunlop et al.l 120031 ) . In order to check the va- 
lidity of our assumption on the host-galaxy, we 
estimate the associated R-band (A'"'^''* = 0.7/iTO) 
absolute magnitudes, M^^, and compare them 
with those measured for the large sample of 

The measured 



quasars m 
Mr range 



Dunlop et al 



(120031 ). 

to -22.2, 



from —25.9 to —22.2, and the me- 
dian value is —23.7. These values a re consistent 
with those found by iDunlop et al. 



= — 23.53±0.09, but our sample shows a wider 
dispersion that can be attributed to the large 
uncertainty of our method and to the different 
redshift range. A host galaxy of later type cannot 
be ruled out, but the low number of detections 
in the optical, the lack of NIR (JHK) data, and 
contamination from the AGN light do not allow 
us to better constrain the host type. In most of 
the cases, the elliptical template provides an ac- 
ceptable fit. Only in one source, LH_A8, the SED 
of the stellar component is not well fit with an el- 
liptical template because of an excess of emission 
at A'''^''* <0.3/im and a late spiral template could 
provide a better fit. Note that a heavily obscured 
starburst galaxy would have an optical-NIR SED 
that is very similar to that of an elliptical galaxy. 

The reduced obtained by comparing the 
model with the data at A > Ifim in the rest-frame 
are reported in Table [H We will refer to these 
as T (for torus) models (see also Table |4|). In 11 
cases, a T model well reproduces the observed IR 
SED and spectrum {x1 <1). In the remaining 10 
cases, poorer fits are obtained (1< xt <10). In 3 
cases the model underestimates the NIR emission 
(LH_A4, LH_A5, and N2_08), in 4 cases the model 
fails to reproduce the depth of the Si feature (Bl, 
B4, MPS78, and MIPS22303), and in the remain- 
ing 3 cases the large xl values are caused by ei- 
ther a noisy IRS spectrum (LH_02, and N2_09), 
or a poor fit at short wavelengths (N2_09, and 
MIPS22304). There are also 2 sources with good 
xi but their fits do not well reproduce the Si fea- 
ture. This is due to the fact that the feature is 
not well sampled by the representative values for 
the IRS spectrum. There is indeed only one data 
point in the Si feature. 

The main cause of failure in the fits with large 
xl is the inability of the models to reproduce 
both a prominent NIR emission and a deep Si 
absorption feature. Yet, sources with deep Si fea- 
tures are expected to have weak NIR emission (e.g. 
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Pier fc Kroli"^ll992HLevensoii et al.ll2006l ). A NIR 

excess requires that hot dust is seen directly by the 
observer (typicaUy happening only in less obscured 
objects), while the Si absorption feature requires 
the presence of a significant amount of optically 
and g eometrica lly thick dust (e.g. Pier fc KrolikI 
19921: iGranato et all Il997t iLevenson et al.i l2006r 



alJl 

i), 



Imanishi et al. 20071). and this cooler dust com- 



ponent should absorb the hot dust emission from 
the vicinity of the obscured AGN. 

We attempted, therefore, to model the objects 
with poor fits, as well as the rest of the sample 
with an alternative model that includes a cold 
absorber detached from the torus (T+C models 
hereinafter). We assume the Ga lactic Center ex- 
tinction curve ( Chiar et al. 2006f) for the cold ab- 
sorber. We do not model re-emission from this ab- 
sorber for simplicity and because this is expected 
to occur at FIR wavelengths. The number of de- 
grees of freedom in the xt estimates is 4 in case of 
T models and 5 in case of T-l-C models. The cold 
absorber is able to produce deeper Si features, 
and significantly improves the fits for 9 sources, 
as shown in Figure [3l In 2 cases, for B2 and B5, 
the xt have actually increased, but they are still 
<1. In these two sources the xt rnore sensitive 
to the broad-band photometric data than to the 
IRS spectrum because of our limited sampling in 
the observed deep Si feature. But the two fits are 
equally good and from a visual inspection of the 
fits, we choose the model with the cold absorber 
component as preferred model because it better 
fits the Si feature. 

The T-fC models were adopted for these 2 
sources and for 7 sources for which the T-fC model 
gives a significantly better fit (LH_A4, LII_A5, 
N2.08, Bl, B2, B4, B5, M1PS78 and M1PS22303). 
In two cases (B4, and N2_08), even the T-l-C model 
does not provide a good fit to the NIR emission, 
although it is an improvement with respect to 
the T model. No improvement was obtained for 
the 3 sources with poor xt (LH_02, N2_09, and 
MIPS22204), therefore for these we keep the sim- 
ple T model. In addition to the T-l-C models, we 
also show in Figure [3] the torus model before ap- 
plying the extinction due to the cold absorber (see 
dot-dashed curve in Figure [S]). The cold absorber 
can absorb up to 84% of torus emission in the near- 
and mid-IR (1-50 /im), as shown in Figure [31 

In modeling the SEDs, we neglected the optical 



data. However, we notice that in most of the cases 
the predicted optical flux from the torus model is 
lower than the observed optical emission. The op- 
tical light might be associated with the host galaxy 
or with scattered light from the AGN. In 6 sources 
detected in the optical, the SED shows an upturn 
or a blue continuum toward shorter wavelengths, 
i.e. in the rest-frame far-ultraviolet. This feature 
can be reproduced by either a population of young 
stars or by emission from the AGN accretion disk. 
The latter scenario is favored by the optical spec- 
trum of source LII_A4 (aka SW104409). Its opti- 
cal spectrum is dominated by a blue faint contin- 
uum, narrow emission lines, e.g. Lya and C IV 
A1549 with asymmetric an d weak broad compo- 



nents ( Polletta et al. 20061 ). The most likely ex- 



planation for the properties of the observed opti- 
cal spectrum is that it is dominated by scattered 
light, and the scatterin g frac tion is estimated to 
be <1% ()Polletta et al.ll2006h . 

In one case (N2_06), the LOS is at only 30° from 
the torus axis, implying an almost clear view of the 
nuclear region. Since no additional absorber is re- 
quired by the data, we expect to see the nuclear 
emission directly. This seems in apparent con- 
tradiction with the faint optical emission of this 
object. Its best-fit model, however, predicts one 
optically thick cloud, associated with the torus, 
along the LOS (see Table [4|). The presence of 
such a cloud in front of the optical source would 
be enough to absorb the intrinsic optical emission 
in this source. The estimated apparent Si optical 
depth in the model is TSi=0.07-0.34 which corre- 
sponds to a visual optical depth of 0.4-2.1 or to 
a suppression factor of the optical flux of 40%- 
0.8%. However, Tsi is poorly constrained by the 
data and in the model because this source is at 
z=2.96. The red MIR SED of this object would 
imply a much higher optical depth than derived by 
the Si feature. Because of its high z, the dust dis- 
tribution in N2_06 is thus poorly constrained. A 
deeper Si might be present implying a higher torus 
inclination or the presence of a cold absorber. On 
the other hand, it is also possible that our best-flt 
model is correct and that the source has only a 
weak Si absorption feature and a low inclination 
torus, and one optically thick cloud is responsi- 
ble for suppressing the optical and NIR nuclear 
light. This scenario might be quite common in 
obscured QSOs as suggested by a recent study of 
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Fig. 3. — Rest-frame spectral energy distributions in uLi, vs A (black symbols and curve) and best-fit model 
(solid grey curve) obtained from a torus-|-cold absorber model (dashed grey curve) and an elliptical template 
(dotted grey curve) to fit the residuals in the near-lR. The dot-dashed curve represents the torus model, in 
the 1-50/xm rest-frame wavelength range, before including the cold absorber. Symbols as in Figure [H The 
name and redshift of each source are annotated. 



the IR S spectra of a sam ple of type 2 absorbed 
QSOs (ISturm et al.ll2006l) . 



5.1. Model parameters 

Here, we analyze the model parameters of the 
best-fit models. In case of the 9 sources for which 
the T-f-C model is preferred, we consider only the 
parameters obtained with such a model (see Ta- 
ble S]). 

The cloud density distribution, which can be 
considered as an indicator for a compact or ex- 



tended torus, is expressed by a power-law, nr{r) oc 
r^°, where r is the torus radius. In the models 
grid, the index a varies from 1 to 3 in steps of 
0.5. Larger values of the index a indicate more 
compact distributions. Since values lower than 1 
are not supported by theoretical consideration s for 
an accretion scenario iBeckert fc Duschlll2004l ). we 
set a minimum value for the parameter a of unity. 
In models with a >1, the depth of the Si feature 
increases for larg er values of a (see Fig. 10 in 



H5nig et al.ll2006l i 



Trends of Si strength with a are different for 
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our models based on an accretion torus co mpared 



to the models in Levenson et al 



They 

do not model the dust distribution in a torus but 
instead use different dust distributions (slab and 
shell) around a central heating source and they do 
not consider the combination of a clumpy medium 
and an ad ditional absorbe r as pr oposed here. Ac- 
cording to iLevenson et alJ (j2007l ). the strength of 
the Si absorption feature is a function of the tem- 
perature contrast in dust material, and, thus, a 
deep absorption feature favors a more extended 
dust distribution, i.e. a=0-l. However, for the 
majority of our sources, a more compact dust dis- 
tribution is favored; a=3 is chosen in 9 cases, a=2 
in 10 cases, and smaller values, a <1.5, are pre- 
ferred in only 2 sources. 

The clouds' vertical distribution (flaring or non- 
flaring) is approximated by a power-law, H(r) oc 
r'', where H is the torus scale height and r the 
torus radius. In the models grid, the index b can 
assume values 1, 1.5, and 2. A non-flaring dis- 
tribution (6=1) is preferred by 17 sources, and 
moderately flaring (6=1.5) is preferred only by 4 
sources. 

The torus inclination is defined by the angle 9 
between the torus axis and the LOS. A face-on 
torus has 6—0° and an edge-on torus has 6=90°. 
In the models grid, the angle 6 varies from to 
90° in steps of 15°. Our sources do not show a 
preferred torus inclination, indeed 6 ranges from 
15 to 90°, with most of the sources (13 sources) at 
intermediate values (6'=30"45°). AU but two (B5, 
and M1PS22303) of the sources modeled with the 
T-l-C model and one of those with weak Si ab- 
sorption feature (N2_06) favor a torus with little 
inclination, 6'=0-30°. All of the others have in- 
clined tori, 6 >45°. Therefore, optical obscuration 
or Si in absorption do not neces sarily imply a LOS 
intercepting the t orus (see also Rigbv et al. 2006t 



Brand et al.ll2007 ) 



An additional model parameter is the number 
of clouds along the LOS, N^'^^ . This parame- 
ter is the average number of clouds obtained from 
five different cloud arrangements for the same set 
of model parameters, and is thus indicative of the 
extinction to the A GN emission produced by the 



torus, ~ TVg^'^^ (|Natta fc Panagialll984f) . The 



apparent optical depth in the Si feature of the T 
and T-t-C models, t|J, and are measured 

following the same procedure applied to the data 



(see § [4]) . The apparent optical depths so derived 
are reported in Tabled) The line-of-sight optical 
depths associated with the torus in the T models 
are TgJ <1.04, with a median value of TgJ =0.13 
or 0.46 depending on the method applied to esti- 
mate TgJ. The optical depths associated with the 
cold absorber are ^^=4-25, with a median value 
of 15. The apparent optical depths for the sources 
fit with the T+C model, t^^'^ , range from 0.6 to 
2.9, with a median value of 1.6, or from 0.4 to 
1.4 with a median value of 0.9 depending on the 
method used to measure it. 
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Fig. 4. — Comparison between the apparent Si op- 
tical depth at 9.7/im measured from the SED and 
spectrum, as described in § [2l and from the torus 
models, as described in § [5l Full circles represent 
sources fit with a T model, diamonds represent 
sources fit with a T-l-C model. The dashed line 
represents the best linear fit between the modeled 
and the measured Tsi, and the dotted lines include 
the la uncertainties on the best-fit parameters. 

In Figure [H we compare the modeled optical 
depths in the 5*1 feature, r^f^, with the observed 
optical depths, r'^gf"-" (see § |4]). The modeled 
optical depths range from 0.01 to 2.9, consider- 
ing both the T and T-f C models. The modeled 
and measured rsi satisfy the following relation, 
T™°'^- =0.96±0.13xT™'=""- - 0.03±0.13. Although 
there is a general agreement between the two esti- 
mates, it is clear that these kind of measurements 
are characterized by large uncertainties, especially 
in objects at high redshift for which the Si feature 
is not well sampled, and that it is a difficult task 
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to model both the NIR-MIR continuum and the 
Si absorption feature. 

In summary, we find that there is no preferred 
torus inchnation associated with the detection of 
a Si feature in absorption. However, in the cases 
where the absorption feature is well detected and 
can be modeled by the torus model, the inclination 
is always higher than 6—45°, consistent with the 
LOS intercepting the torus. Sources with a less 
inclined torus require an additional absorber to 
explain the observed Si feature or show a weak Si 
feature. A compact non-flaring torus is preferred 
by the majority of the sources. The clouds' ra- 
dial density distribution indeed indicates that the 
torus emission region is compact, or that the near- 
and mid-IR emission is dominated by dust in the 
vicinity of the nucleus. The preferred non- flaring 
torus in the majority of these luminous sources 
is consistent with the predictions of the receding 



IS consistent witn tne predictions oi tne recedmg 
torus models (e.g. ^Simpsonl l2005 : Honig fc BeckertI 
20071 ). According to these models, the opening 



angle of the torus increases at larger luminosi- 
ties. Our result s are in ag reement with the pre- 
dictions from , Honig fc Beck ert (2007). These au- 
thors claim that flaring should not occur in high 
luminosity sources because large clouds at large 
distances from the AGN should be driven away by 
the radiation pressure. 

5.2. Far-IR emission 

According to current AGN evolutionary mod- 
els, obscured and extremely luminous AGNs are 
believed to represent a specific and rare phase in 
the evolution of an AGN, when the central super- 
massive black hole (SMBH) is fully grown and 
still surroun ded by a larg:e amount of dust and 
gas (San ders et all Il988t IPi Matteo et all l2005t 
iHopkins et al.ll2005f l. High luminosity AGNs are 
believed to be triggered by large scale galaxy 
mergers and, therefore, to be accompanied by in- 
tense starburst activity. In order to test whether 
this scenario applies to our sample, we search for 
starburst signatures in our objects. Obscured star- 
burst galaxies are generally heavily extincted at 
optical wavelengths and are weak X-ray sources 
compared to AGN. Thus, the only observations 
that could reveal a starburst are PAH features at 
MIR wavelengths and strong continuum from cool 
dust at far-lR (FIR) and sub- mm wavelengths. An 
analysis of the starburst contribution based on the 



PAH features is discussed in section 19.11 here we 
analyze the FIR properties. 

Most of the power produced by a starburst 
emerges in the FIR. The typical FIR emission 
of powerful starbursts is characterized by lumi- 



nosities ~10 



11. 0-12. 5t 



and peaks at 60-100 



/^m (jSanders fc Mirabell 1996!). Such luminosi- 
ties correspond to star- formation rates (SFRs) of 



~20 to 550 Mq yr-^ (|Kennicutd [l99al . Higher 
rates, up to ^^SOOOM© yr~^, and luminosities, 
up to 3xlO^'^L0 are measured in some high- 2: 
ULIRGs/sub-mm galaxies (SMGs) and a ttributed 
to starbursts (e.g. iGhapman et al.ll2004l ). 

In our sample, 5 SWIRE sources are detected 
at 70 or 160/xm and 3 E-FLS sources are detected 
at 70^m. Note that both MIPS FIR observations 
are available only for SWIRE sources (13 sources), 
and MIPS 70/im pointed o bservations are avai l- 
able for the E-FLS sources ( Sajina et al.ll20()7ah . 
The SWIRE observations in the FIR are sensitive 
enough only to reveal FIR luminosities >10^'^Lq 
at the observed redshifts (the 5a limits are 18 mJy, 
and 108 mJy at 70, and 160/xm, respectively). In 
Figure El we compare the SWIRE 5a detection 
limit and the detected fluxes at 70 and 160/im of 
our sources with those expected for two starburst 
galaxies, M82 and Arp220, with an IR luminosity 
of 10"Lo and IO^^.Sl^ ^ ^3 2. The figure 

clearly shows that we do not expect to detect the 
FIR emission of starburst with such luminosities 
at z >1.2. 

What is then the origin of the detected FIR 
emission in our objects ? We investigate three 
possible origins, the torus, AGN-heated dust at 
large distances from the nucleus, and an excep- 
tionally powerful starburst. The spectrum of our 
torus models peaks at around 10/zm, correspond- 
ing to a temperature of about 300 K, and falls off 
at longer wavelengths. In two cases, the 70/im de- 
tections or upper limits agree well with our models 
(Nl_09, N2_06), but in the other six cases, the ob- 
served 70 or 160/im fluxes are significantly higher 
than the model predictions (LH_A8, N2_08, N2_09, 
M1PS42, M1PS22204, and MIPS22303) and re- 
quire an additional component to be explained. 
The SEDs of these sources resemble that of the 
starburst/ AGN composite s ource CXO-J1417 dis- 
cussed in iLe Floc'h et"al] (|2007t ). The FIR lu- 
minosity and SFR estimated for CXO-J1417 are 
4.5xlO^^L0, and ~750MQyr~^, respectively, as 
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Fig. 5. — Predicted MIPS 70/j,m (panel a) and 160/im (panel b) fluxes for the starburst galaxies M 82 and 
Arp 220 with IR luminosities equal to IO^^Lq (dotted curves) and lO^^'^L© (solid curves). The horizontal 
dashed line shows the SWIRE 5a limit. Full black circles and diamonds represent sources in the selected 
sample detected at 70/im or 160/im fit with a T and a T+C model, respectively. Source names are annotated. 



in a powerful starburst. 

In order to quantify the amount of luminos- 
ity detected in the FIR, we model the observed 
fluxes at rest- frame A >5/xm with a starburst tem- 
plate. We first derive the residual fluxes after 
subtracting the torus contribution to the observed 
FIR fluxes and then fit them with starburst tem- 
plates. Two different starburs t templates, M82 

1998h . are used for the 
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and Arp 220 

fits, but we only show the results obtained with 
the M82 template because it provides better fits 
in the majority of the cases and similar good fits in 
the remaining cases. The starburst fits, combined 
with the torus and host galaxy fits, are shown in 
Figure [6l The estimated starburst luminosities 
and contribution to the total IR luminosity are 
reported in Table [31 In half of the FIR-detected 
sources the torus contribution to the total IR lu- 
minosity is lower than that of the additional FIR 
component. The remaining objects divide equally 
between those for which the two contributions are 
similar and those for which the torus contribution 
is larger. The estimated starburst FIR luminosi- 
ties are all greater than 3.3x10^^Lq. These lu- 
mino sities imply SFRs > 600 yr~^ ( KenniciifrO 
1998I ). which are consistent only with those mea- 
sured in the most powerful starbursts. 

Another possible origin of the FIR component 
is re-emission of the energy absorbed by the cold 
absorber discussed in iJSl There are two sources 



detected at long wavelengths and modeled with 
the T-j-C model, N2.08 and MIPS22303. The 
estimated absorbed luminosities in the 1-50 /im 
wavelength range are 1.8x10"*^ ergss"^, and 
5.25 xlO^'s ergs s-\ for N2.08 and MIPS22303, 
respectively. These correspond to 2.8 and 0.24 
times the FIR luminosities derived from the star- 
burst template (iff^). Although these ratios 
depend on the cold absorber covering factor and 
on the efficiency in re-emitting the absorbed en- 
ergy, these values indicate that this hypothesis is 
energetically viable. 

Although we can reproduce the FIR compo- 
nent with a starburst template or with re-emission 
from the cold absorber, its origin remains undeter- 
mined. Large FIR luminosities, as measured in our 
sources, even when cold molecular gas is detected, 
do not probe the presence of a powerful starburst, 
e.g. in Q SOs wher e [O II] A3727 emission is not 
detected (|Holl2005h . Thus, both explanations re- 
main plausible and do not exclude others. In order 
to understand the origin of the FIR component 
we would need more FIR measurements, higher 
signal-to-noise IR spectra, and near-IR spectra. 
Multiple FIR detections would help to constrain 
the radiation field, high S/N IR spectra would al- 
low us to search for and, eventually, measure PAH 
emission (see H9.1\i . and near-IR spectra could be 
used to constrain the SFR from the [O II] A3727 
emission line. 
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Fig. 6. — Rest-frame spectral energy distributions in vLn vs A (black symbols and curve) of sources detected 
at 70 or 160/im. Also shown are the best-fit model (solid black curve) obtained from a torus (-|-cold absorber) 
(dashed grey curve), an elliptical template to fit the residuals in the near-IR (dotted grey curve), and a 
starburst template (dot-dashed grey curve) to fit the residuals in the far-lR (triangles). Symbols as in 
Figure [T] The name and redshift of each source are annotated. 



5.3. AGN bolometric luminosity 

Assuming the best-fit torus model and inte- 
grating it from lOOOA to lOOO/im, we can derive 
a lower limit to the AGN bolometric luminosity. 
This estimate of AGN bolometric luminosity does 
not include the AGN optical light that is not re- 
processed by the dust and it is not emitted along 
our LOS. Thus, the true AGN bolometric luminos- 
ity should be corrected by a factor that depends 
on the torus covering factor and anisotropy of the 
torus re-emission. The covering factor of a torus 
can vary from 20 to 50%, implying AGN bolomet- 
ric luminosities that can be 2-5 times larger than 
our estimates. The uncertainty associated with 
the redshift estimates (±0.2) corresponds to a 30% 
(20%) uncertainty in the luminosity estimate for 
objects at z ~2 (z ~3). We do not include the 
luminosity associated with the host galaxy. Note 
that the discrepancy between the observed optical 
data and the best-fit models does not affect the 
results, since the bulk of the luminosity is emitted 
around lO/Ltm in the rest-frame, where the SED 



is well reproduced by the models. The estimated 
bolometric luminosities, reported in Table[3l range 
from 10"*^ to lO^'^ergss^^ 

The estimated AGN bolometric luminosities are 
compared with the 6/xm luminosities in Figure [7] 
in order to investigate whether it is possible to 
determine AGN bolometric luminosities from the 
observed 6/im luminosities without modeling the 
SED and without having measurements at long 
wavelengths. We find that the two luminosities are 
highly correlated and that the AGN bolometric lu- 
minosity can be derived from the 6/im luminosity 
adopting the following relationship: 

Log{Lf^^) = Log{L{Q^im)) + 0.32 ± 0.06. (1) 

This relationship gives a bolometric luminosity 
that is seven times smaller than what would be 
derived assuming the bolometric corrections ob- 
tained for a sample of o ptical ly-selected, type 
1 QSOs by lElvis et al.l (Il994l ). Log{Lboi) = 
Log(L{6^m)) + 1.16. Such a difference is ma inly 
due to the fact that Lbo; in I Elvis et al. m- 
cludes both the AGN optical emission that is not 
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included in our analysis, and the emission from 
the host galaxy. Their Lboi also include a signifi- 
cant contribution in the FIR from cool dust which 
is not included in our models and thus in our esti- 
mates of L-^^^ . Moreover, their QSOs are about 
100 times less MIR luminous than our sample, 
and, thus, it is also probable that their median 
template does not well represent the MIR-FIR 
SEDs of our objects. 




.6 45.8 46.0 46.2 46.4 45.6 
Log(L(6Mm)) (Lg) 

Fig. 7. — Comparison between the measured 6/im 
luminosity and the estimated AGN bolometric lu- 
minosity derived by integrating the best fit torus 
model between 0.1 and 1000/im. The dashed line 
corresponds to a linear fit with slope 1, and the 
dotted lines correspond to the la dispersion. The 
bolometric luminosities derived by including the 
FIR luminosity in sources with FIR detections cor- 
respond to the top of the solid vertical lines. The 
additional FIR luminosity is derived by fitting the 
residuals in the FIR with an M82 template (see 
Figure [6]). 

For the 8 objects for which FIR data are avail- 
able (see previous section), we also estimate a 
bolometric luminosity including the FIR compo- 
nent (see values in Figure [7]) . In those cases the 
bolometric luminosities are on average 2.5 times 
larger than estimated from the torus model. This 
implies a significant amount of cool dust which 
dominates the bolometric luminosity and which is 
located further from the nucleus than the torus. 
This dust could be illuminated by radiation from 
the AGN that is not intercepted by the torus, 



which would be, qualitatively, consistent with the 
T-l-C models described above. Overall, by taking 
into account the torus covering factor, and a FIR 
component, the bolometric luminosities of our ob- 
jects can be up to 7-17 times larger than estimated 
using equation [TJ 

6. Composite IR spectra 

In order to search for weak features in the IRS 
spectra of Figure 2, the individual spectra, after 
being boxcar-smoothed to three times the approx- 
imate resolution of each module, were combined 
to obtain a composite spectrum with higher S/N. 
Composite spectra were made for 2 sub-samples: 

1) sources fit with the T model (12 sources); and 

2) sources fit with the T-f-C model (9 sources). 
The composite spectra were obtained by taking 
the median of all available data in the individual 
spectra in bins of A{Log{X))=0.005 with A in /im, 
after normalizing them at 6fim in the rest-frame. 
The width of each bin is increased if less than 3 
measurements are available and it is reduced if 
more than 100 measurements are available. The 
uncertainty associated with the median value per 
bin corresponds to the uncertainty of the mean. 
The two composite spectra and associated la un- 
certainties are shown in Figure [H The composite 
spectrum of the T group is characterized by a large 
dispersion, especially at A ~10/im. This large dis- 
persion indicates that we combined sources with a 
broad range of depths in the Si feature. 

Some spectral features other than Si that might 
be detected in the two composite spectra are high- 
lighted, e.g., absorption features due to ice at 
6.15/xm and hydrocarbons (HAG) at 6.85^m, the 
atomic emission features [ArII]A6.99/im, [NeVI]A7.64/xm, 
and [SIV]A10.54/im, the molecular hydrogen emis- 
sion line H2(S3)A9.66/xm, and the 6.2, 7.7, 8.6, 
and 11.3fim PAH emission features. Of these, 
the strongest feature that may be detected in 
both composite spectra is the 7.7/im PAH. This 
PA H is the strongest feature in starburs t galax- 
ies (|Sturm et al.ll2000l: iBrandl et al.ll2006l) . 

The strongest feature observed in the compos- 
ite spectra is associated with molecular gas H2, 
H2(S3)A9.66/im. The detection of such a feature 
cannot be confirmed because of the limited signal- 
to-noise, but would likely indicate the presence 
of shocked molecular hydrogen outside the MIR 
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Fig. 8. — Median normalized IRS spectra of two 
sub-group of sources: (a) 12 sources modeled with 
the T model; and (b) 9 sources modeled with T+C 
model. The shaded area represents the mean ab- 
solute deviation. The dashed curve represents a 
power-law of slope a— 2 (F^ oc v~°') normalized at 
6/im. Some spectral features that might be present 
are annotated. 



absorbing material, as can be expected in case 
of high- velocity galaxy collisions, starburst super- 
winds, or AGN-driven outflows. Further observa- 
tions would be necessary to test these scenarios. 

Absorption due to HAC might be present 
in both composite spectra. This might be 
associated with the material that is also re- 
sponsible for the deep Si fea ture, as in IRAS 
00183-7111 (ISpoonet al.ll2004l) . 

In order to facilitate a comparison between the 
two spectra, we overplotted a power-law model 
with slope of 2 normalized at 6/xm in both panels 
of Figure [SI The main differences among the spec- 
tra are the 2~6/im spectral slopes and the depths 
and shapes of the Si feature. In group (1) (T 
model), the FWHM of the Si feature is 1.8/xm, the 
apparent optical depth is 0.4, and the 2-6/xm rest- 
frame power-law slope is 3.0. In group (2) (T-|-C 
model), the FWHM of the Si feature is 3.9/zm, 
the apparent optical depth is 0.9, and the 2-6/im 
rest-frame power-law slope is 1.7. It is interest- 
ing to note that, contrary to the expectation from 
torus models that the steepest NIR spectra are as- 
sociated with the shallowest Si absorption feature, 
the spectrum with the deepest Si feature (group 
2) shows instead the flattest NIR slope. This re- 
sult supports the scenario proposed in § [5] for the 
sources fit with the T+C models in which the hot 
dust close to the AGN is directly visible, and an 
additional absorber is required to reproduce the 
deep Si feature. 

7. Comparison with other AGN samples 
7.1. Fraction of obscured QSOs 

In order to estimate the total surface density of 
the QSOs in our IRS sample, we first define em- 
pirically the colors that encompass our IRS sam- 
ple. This is done by analyzing the location of our 
sources in various IR color-color diagrams. Then, 
we select a reference flux-limited sample of IR 
sources from a large and contiguous area which oc- 
cupy the same location in the selected color-color 
diagrams. To select our reference sample we used 
the SWIRE catalog in the Lockman Hole field. 

It has been established that the majority of 
sources with red IRAC colors from Spitzer sur- 
veys have the colors expected f rom SEDS o f AGN 



at va r ious redshifts (iLacv et al. 2004 : Stern et al 



2005t iHatziminaoglou et al. I2OO5I: iBarmbv et al 
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Fig. 9. — Panel a: IRAC color-color diagram. Panel b: IRAC-MIPS color-color diagram. The sources 
discussed in this paper are shown as red stars. Black circles represent all of the sources in a randomly selected 
ldeg2 field with F (24/im)> 209m Jy and F(8Mm)>47MJy- The gr een fuU squares represent the SDSS/SWIRE 
type 1 QSOs with 1 < z < 3 from iHatziminaoglou et ah ( 2005h . The curves show the color of templates of 
various types as annotated (solid: 0.1< z <2, dashed: 2< z <3). The regions to the right of the dashed 
lines represent the color criteria applied to select sources with colors consistent with those in our sample: 
unobscured QSO candidates (cyan triangles) and obscured QSO candidates (magenta downward triangles) 
(see §[ni). 



20061 ). AGNs have such colors because of their 
hot dust emission. We thus analyzed several color- 
color diagrams combining the 4 IRAC bands and 
the MIPS 24 fj,m band. From this analysis, we 
conclude that we can select sources with similar 
colors as our IRS sample by using only two color- 
color diagrams. One diagram combines all 4 IRAC 
bands and it is shown in the left panel of Fig- 
ure [HI The other diagrams is made by combin- 
ing 3 bands, IRAC[3.6], IRAC[8.0], and MIPS[24], 
and it is shown in the right panel of Figure [9l 
In both panels, we show our IRS sample and 1540 
sources in a 1 deg^ randomly selected field from the 
SWIRE survey. The 1540 sources are all detected 
in 4 IRAC bands and at 24Mm to the SWIRE lim- 
its in the Lockman Hole field (4.2, 7.5, 46.1, 47.3, 
and 209fi3y at 3.6, 4.5, 5.8, 8.0 and 24:fim, respec- 
tively). 

The highly luminous and obscured QSOs in 
our IRS sample occupy the areas in the color- 
color diagrams with the reddest colors. The 
only template that shows similar colors corre 



spon ds to a highly obscured AGN (jPolletta et al 
I2OO6I ). Based on the areas in the diagrams 



occupied by our IRS sample, we define the 
following multi-color criteria to locate similar 
sources: F(5.8)/F(3.6)>2, F(8.0)/F(4.5)>2, and 
Log(F(8.0)/F(3.6))-l-Log(F(24)/F(3.6)>2.0. These 
color criteria are also shown in Figure [9] by dashed 
lines. 

Having defined this color space for the highly 
luminous, obscured QSOs in our spectroscopic 
sample, we select from a random 3.5 deg^ field 
within the SWIRE Lockman Hole field all sources 
detected in the 4 IRAC bands and with 24/im 
fluxes greater than 1 mJy with colors consistent 
with the color criteria mentioned above and shown 
in Figure [51 

Using all of these criteria, we find a total of 
82 sources in 3.5 deg^ in this random sample of 
SWIRE sources. To evaluate in more detail if all 
of these sources are similar to our obscured QSOs, 
we first visually inspected all of the available im- 
ages to remove any source with spurious colors. 
We found 1 confused source, 1 saturated star and 
2 nearby galaxies with bad aperture photometric 
data. Thus the sample was reduced from 82 to 78 
sources (22 deg~^). Next, we examined the over- 
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all SEDs of the 78 sources, including the optical 
data. Based on this analysis, we estimate that 
all selected sources have SEDs consistent with be- 
ing AGN dominated and with those in our sam- 
ple. We have previously found that sources clas- 
sified as starbursts by IRS spectra can be distin- 
guished from AGN in the SWIRE sam ples based 
on their SEDs ( Weedman et al.l[2006al ). so we can 
have confidence that the SED analysis does cor- 
rectly locate AGNs. 

To further test that our criteria select only 
AGN-dominated sources at z=l. 3-3.0, we ap- 
pli ed them t o the large sample of 48 sources 
in lYan et all (|2007l ). for which a MIR classifi- 
cation, Spitzer IRAC and MIPS 24/im data, and 
spectroscopic redshifts are available. Our crite- 
ria select 7 sources from this sample, of which 
4 are included in our sample. The remaining 3 
(MIPS8034, MIPS15958, and MIPS22277) are all 
AGN-dominated in the MIR with EWy.T^™ <0.22, 
two fall into the rcdshift and luminosity range of 
our sample, but one (MIPS8034) is at lower red- 
shift (2=0. 95) and lower luminosity (L(5.8/im) = 
1.6x10"'^^ L©)- Note that the other two sources are 
not formally included in our IRS sample because 
MIPS15958 has a MIR luminosity of IO^-^^Lq, 
that is just below our selection threshold, and the 
IRS spectrum of MIPS22277 shows some PAH 
emission. This test indicates that our sample is 
highly reliable in selecting AGN, but it might in- 
clude some (1/7 or 14%) sources at lower redshift 
and MIR luminosities. 

We can perform an additional test on the de- 
gree of contamination from low redshift and low 
luminosity sources using the available optical spec- 
troscopic data in the selected sample. Among the 
78 selected sources, 12 sources have optical spec- 
troscopic data. All 12 sources with spectroscopic 
data have redshifts from 1.3 to 3.0, with a median 
value of 2.5. This is the same range of redshifts 
found for our IRS sample, indicating that our se- 
lection succeeds in selecting sources in a similar 
redshift range. If 1 out of 7 sources in our sam- 
ple, or 14%, were at lower 0, a s suggested by the 
analysis of the lYan et al.l (j2007l ) sample, we would 
have expected 1-2 low-z sources in our spectro- 
scopic sample. Since no low-z, low-L interlopers 
were found in the spectroscopic sample, we can 
set a limit on the contamination fraction of <8%, 
which corresponds to 1 out of 12 sources, the min- 



imum number we could have found in our spectro- 
scopic sample. A redshift greater than 1.3 and the 
bright 2Afim flux assure that also the MIR lumi- 
nosity is high and consistent with those measured 
in our sample. 

Since our selection of the 78 sources is only 
based on IR colors, it is possible that some of these 
sources are not as obscured as our IRS sample. We 
use the F(3.6/im)/F(optical) flux ratio, where for 
the optical we used the three bands g', r', and 
i' , to separate obscured and unobsc ured candi- 
dates. Based on previous studies 
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20061 ). we consider obscured candidates those with 
F(3.6/xm)/F((7', r' , and i') greater than 15, 13, and 
10, respectively, and unobscured candidates those 
with at least one flux ratio lower than those thresh- 
olds. This criterion yields 61 obscured AGN and 
17 unobscured AGN candidates. The optical- IR 
flux ratios of the QSO candidates are compared 
with those of our IRS sample and with other sam- 
ples from the literature in Figure [TOl 
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Fig. 10.— IR/optical flux ratios of the MIR- 
selected QSO sample (red stars) and of the 
SWIRE sources with similar IR colors (triangles) 
(selected as shown in Fig.[S]). The F(3.6/im)/F(r') 
flux ratio shows one of the criteria applied to sepa- 
rate obscured (magenta downward triangles) and 
unobscured (blue upward triangles) QSO candi- 
dates. Symbols as in Figure [9l 

In order to verify our simple classiflcation, we 
examined the optical spectral properties of the 
sources with available optical spectra and spec- 
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troscopic classification. A spectral classification is 
available for six sources, and they are all AGN ac- 
cording to the observed emission lines; 2 show nar- 
row emission lines and 4 show broad emission lines. 
Our classification agrees with the spectroscopic 
classification, in the sense that those with nar- 
row lines are classified as obscured and those with 
broad lines as unobscured, in 5 out of 6 cases. The 
only exception is a source with Lya, N V A1240, 
and C IV A1549 emission lines with FWHM-1700 
km s""'^, intermediate between a narrow and a 
broad emission line AGN. Although this compar- 
ison confirms our classification, we want to stress 
that our method based on IR/optical flux ratios is 
based only on a handful of sources whose availabil- 
ity of optical spectra is biased in favor of AGNs. 
Thus, any selection based on the criteria men- 
tioned above should be applied with such a caveat 
in mind. 

Based on this result alone, we would conclude 
that the surface density of obscured QSOs to the 
limits of the SWIRE survey is 17 deg^^ (61 in 
3.5 deg^). However, there are a few sources with 
low F(3.6/j,m)/F(optical) flux ratios, but with ex- 
tremely red IR SEDs where the 3.6/im flux is much 
lower than the fluxes at longer wavelengths. Thus, 
it is likely that some unobscured QSO candidates 
might actually be obscured or reddened type 1 
AGNs. In such a case, we might be underesti- 
mating the surface density of obscured QSOs. 

Indeed our color selection criteria rule out most 



of the type 1 QSOs such as those in Hatzimin aoglou et al 
{^200^ . The IR colors for the SDSS/SWIRE type 
1 QSOs in the 1.3-3.0 redshift range are also 
shown for comparison in Figure [9l While there 
is a large overlap between the SDSS/SWIRE 
type 1 QSO sample and our obscured sample 
in IRAC colors (left panel in Figure [9]), the 
SDSS/SWIRE type 1 QSOs are all clustered to- 
wards the bluest side of the color distribution. 
Moreover, there is almost a complete separation in 
the F(24)/F(3.6) color between the two samples, 
with all the SDSS/SWIRE type 1 QSOs having 
F(24)/F(3.6)<30, and all but one of the obscured 
QSOs having F(24)/F(3.6)>30. This difference 
is likely due to the large NIR extinction affecting 
our sample which is absent in the SDSS/SWIRE 
type 1 QSO sample and highlights the ability of 
this flux ratio to distinguish obscured from unob- 
scured QSOs. Thus, our color selection criteria 



disfavor the classical type 1 QSOs. If we consider 
all selected sources as obscured AGN candidates, 
their surface density increases to 22 deg~^ (78 in 
3.5 deg2). 

In order to derive the fraction of obscured 
QSOs among the 24/im bright AGN population at 
z=l. 3-3.0 and with MIR luminosities >1O^2L0, 
we need to estimate the surface density of unob- 
scured, type 1 QSOs with the same redshift range 
and MIR luminosities as our IRS sample. We es- 
timate such a surface densit y from the 8^m lu- 
minosity function derived bv iBrown et al. ( 20061 ) 
for aU AGN with F(24/im)> lmJy. Assum i ng th e 
same QSO template as in I Brown et al. ( 20061 ). 
the derived 8/im absolute magnitude correspond- 
ing to L(6Aim)=lOi2L0, is Mg^™ = -28.39. The 
estimated surface density of all type 1 AGNs 
with Mgpm < —28.39 in the redshift range 1.3< 
z <3.0 is 11.7 deg"'^. A value only slightly 
smaller, 9.4 deg^^, is derived assuming the op- 
tical lu minosity function for t ype 1 AGNs mea- 
sured bv lRichards et all (|2006l ) and Mb < -22.3, 
which corresponds to Mg^^ < —28.39 assuming 
the same template. Using the result derived from 
the MIR-selected QSO sample, the ratio of ob- 
scured to unobscured MIR luminous QSOs in the 
1.3-3.0 redshift range is, therefore, 17-22/11.7, or 
~1.5-1.9. This value is consistent with the frac- 
tion of obscured AGNs measu red in Mill-radio se- 
lected samples (Martfnez-Sansi gre et al. 1|200.5|). 

7.2. Comparison with other IRS MIR- 
selected samples 

As described in § [21 the main selection cri- 
teria applied to carry out IRS observations of 
Spitzer sources were bright 24/xm fluxes and faint 
optical magnitudes ( Houck et al.| [2005). In addi- 
tion, other selection criteria were applied in some 
cases to favor the selection of starburst or AGN 
systems. In order to favor starburst-dominated 
galaxies, lYa^ et al 



( 20061 ) als o require large 



F(24/xm ) /F(8/xm) flux ratios. IWeedman et al 
(|2006ah require an X-ray detection to select AGN 
and a bump in the IRAC SED to select starbursts. 

In FigurefTTl we compare the F(24/im) /F(8/im) 
and F(24/im)/F(r') flux ratios of the sources in 
our IRS sample with the selection criteria ap- 
plied to select the sample from w hich a large 
fraction of our sources was drawn (Houck et al 
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Fig. 11. — Color-color diagram showing the 
selection of the MIR-selected samples from which 
some of the sources in this work are drawn. 
The dotted Hne shows F(24A(m)>0.9mJ y and 
mag(r')>22 from IWeedman et all ( 2006al ) (note 
that they also require an X-ray detection). 
The dot-dashed line sho ws F(24Mm)> 0.75mJy 
and mag(R)>24.5 from iHouck et al.1 (j2005). 
The region enclosed by the dashed lines 
shows Lo^(iyF(2AiJ,m)/iyF(8iJ.m))> 0.5 an d 



Lo£(j/F(24/xm)/i/F(r'))>1.0 from lYan et al 
( 2007t ). For comparison we show all the 
SWIRE sources with F(24^m)>lmJy in 1 
deg2 (open black circles), the SDSS/SWIRE type 
1 QSO from [Hatziminaoglou et al.1 (l2005l) fgreen 
square s), the MIR-selected sources in lYan et al. 
(|2007t) (brown diamonds), and the SWIRE 



sources detected at 8 and 24/xm selected as shown 
in Figure [H (downward magenta triangles and 
blue triangles). Symbols as in Figure [51 



20051: 1 Yan et"aIll2006HWeedman et al.ll2006al) . We 

al so show for compa rison the colors of the sources 



Yan et aP ( 20061 ). of those selected by the IR 



criteria described in the previo us section, of the 



SDSS/ SWIRE type 1 AGN from lHatziminaoglou et al 



(l2005l) . and of a random sample of SWIRE sources 
with F(24/^m)>lmJy in 1 dcg^. 

It is interesting to note the separation in 
F(24/im)/F(8/xm) flux ratios between the selected 
samp le and the starburst and composite sample 
from lYan et"aLl (|2006l ). Lower F(24^m)/F(8^m) 
flux ratios favor AGN-dominated sources, but 
high ratios do not exclusively select starburst 
galaxies. Indeed half of our sources show large 
F( 24/im)/F(8 / im) a nd a large fraction of sources 
in I Yan et all (12006) contain an important AGN 
component ( Sajina et al. 2007a[ ). 

Our sample shows the largest F(24^m)/F(r') 
ratios of all samples. The sources selected in the 
previous section to derive the surface density of 
our IRS sample show F(24/im)/F(r') flux ratios in 
between the SDSS/SWIRE type 1 QSOs and the 
IRS sources. This indicates that the sources pre- 
sented here might actually be more extreme than 
the sample selected above. 

7.3. Luminosity, Tsi and redshift 

The MIR (6/im) luminosities and redshifts of 
the AGN sample presented here are compared 
with those of other AGN samples from the litera- 
ture in Figure [T2l The literat ure samples includ e 



X-ray selected type 2 QSOs (|Sturm et al.ll2006h . 



optically and IR selected (SDSS/SWIRE) type 
1 ^SOs with 1 < z < 3 (jHatziminaoglou et al 



20051 ). X-ray selected, optically faint (r >22) and 



IR bri ght (F(24 ^m)>l mJy) sources (jWeedman et al 



2006ar) ■ MIR selected AGNs and composite sources (|Yan et al 
20061 ). and various types of local AG Ns, type 1 



QSO s, Seyfert 1, Seyfert 2, and ULIRGs (iHao et al 
20071 ). Note that some of the sources in our sam- 



ple are in common with those in IWeedman et al 



2006af) (5 sources), and with those in lYan et al 
20071 ) (5 sources). 



Our sources are among the most luminous MIR 
sources in this compilation. The only sources with 
similar MIR lumin osities are the SDSS/SWIRE 



type 1 QSOs from iHatziminaoglou et aL 



and th e most luminous type 1 QS Os from 



Hao et al 



(2007). The composite sources in fYan et a" ( 2007^ 
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Fig. 12. — Panel a: Mid-IR rest-frame) luminosities vs redshift for the obscur ed AGNs in this w ork (red 
full stars), and various AGN samples, X-ray selected type 2 QSOs (red triangles; Sturm et al. 2006h . X-ray 
selected optically faint AGNs (rever sed magenta tria ngles; IWeedman et al. 2006al ). IR-selected AGNs and 
composite sources (brown d iamonds; Yan et al"]l2007t ). SPSS /S WIRE selected ty pe 1 QSOs (green squares; 
Hatziminaoglou et al. 2005 ). and type 1 QSOs (full blue circles: iHao et al.ll2007l) . Panel b: The same data 
shown in the left panel are shown here with, in additio n, the value s for th e Seyfert 1 (green squares), Seyfert 
2 (cyan crosses) and ULIRGs (open black circles) from lHao et al.l (j2007t ). 



and the type 2 QSOs from lSturm et al.l ( 20061 ) are 
on average 3 and 250 times fainter, respectively. 

The MIR (5.5/im) luminosities and the mea- 
sured apparent Si optical depth, Tsi, of our AGN 
sample are compared with those of other AGN 
samples from some of the samples in the above 
compilation in Figure I1 3[ The literature sam- 



ples include ah AGNs in iHao et all (|2007l) : type 
1 QSOs, ULIRGs, Seyfert 1, and Seyfert 2 galax- 
ies. Our sample covers a new parameter space in 
the MIR luminosity-absorption distribution of all 
known AGNs, as illustrated in Figure [131 The 
selected obscured QSOs show the highest MIR lu- 
minosities (1/(5. 5^m) > IO^^Lq), consistent only 
with the most luminous type 1 QSOs and Seyfert 
1 galaxies, and also show large Si optical depths, 
up to Tsi—3.2, consistent with those observed in 
ULIRGs and in some Seyfert 2 galaxies. 

8. X-ray properties 

X-ray data from Chandra are available for 11 



sources, but only 5 are detected ( Weedman et al 



2006aHPolletta et al.ll2006l . LH.02, LH.A4, LH_A5, 
LH_A6, and LII_A8;). Chandra 5 ksec observa- 
tions are available for the 5 sources in the Bootes 




1 2 
Apparent Tj, 

Fig. 13. — Mid-IR (5.55/im rest-frame) luminosi- 
ties vs the Si apparent optical depth at 9.7/im of 
the obscured AGNs in this work (red full stars) ,^ 
and of various AGN samples from iHao et al 



(|2007D : type 1 QSOs (full blue chcles), ULIRGs 
(empty black circles), Seyfert 1 (green open 
squares), and Seyfert 2 (cyan crosses). Rightward 
pointing arrows represent lower limits to Tsi- 
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field (|Murrav et al.ll2005l) . but none of the selected 
sources is detected. We report an upper limit to 
the broad-band (0.3-8 keV) X-ray flux for these 
sources corresponding to the ir 90% completenes s 



limit (10"" ergscm-^s-i) ([Murray et al.ll2005[ ). 



LH_A11 is part of the 70 ks ec SWIRE / Chandra 



survey ( Polletta et al. 20061 ) . but was not de- 



tected, thus we report an upper limit to the broad 
band flux of 10~^^ ergscm"^ s^^. Chandra ob- 
servations for B2 are scheduled (CXO Proposal 
08700396) and XMM-Newton observations are 
pending for N2_06, Bl, N2.08, and B3 (XMM- 
Newton Proposal 050326). In Table El we hst 
observed X-ray fluxes, estimated effective col- 
umn densities and luminosities for the detected 
5 sources, and upper hmits to the fluxes for the 
remaining 6 sources. Here, we compare the ab- 
sorption measured in the X-rays with the extinc- 
tion measured in the IR and use the MIR to X-ray 
luminosity ratio to investigate the dust covering 
fraction in the selected sample. 

8.1. X-ray absorption and IR extinction 

The effective gas column densities responsible 
for absorption in the X-rays are derived from the 
X -ray hardness r atios using the method described 
in lPolletta et al.l ( 2006 ). The estimated values are 
reported in Table [H The absorption measured in 
the X-rays and the dust extinction measured in the 
IR, Tsi, are compared in FigurefTil We also show a 
curve (dotted curve) corresponding to equivalent 
Nfj and rsi assuming a gas-to-dust ratio of 100, 
instead of the Galactic value, a dust density of 2.5 
g/cm^ density and O.l^m radius per dust grain, 
or Nh^2.26x1022 x Tsi cm-"^. 

In four sources, the Nh is significantly larger 
than that derived in the IR. The only exception, 
where the Nh is lower, is LH_A2 which is charac- 
terized by a soft X-ray spectrum. However, there 
is evidence, based on the X-ray/IR luminosity ra- 
tio, that in LH_A2 the observed X-ray spectrum 
corresponds to a small fraction of the intrinsic X- 
ray emission due to scattering, while the intrinsic 
emission is absorbed by a Compton-thick column 
density (Nh> 10^4 cm-2). This source is also char- 
acterized by the highest IR extinction. 

The large Njj values compared to the extinction 
measured in the IR by rsi are easily explained. 
The apparent Si optical depth underestimates the 
true extinction in the IR because its measurement 
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Fig. 14. — Comparison between the apparent Si 
optical depth and the column density estimated 
from the X-ray hardness ratio. The dotted line 
corresponds to the expected agreement in case 
both absorptions occur in the same medium and 
assuming a Galactic gas-to-dust ratio of 100. The 
dashed curve corresponds to the linear fit found fo r 
a large compilation of AGNs by lShi et~al ( 20061 ). 



does not take into account re-emission of radiation 
at the same wavelengths where absorption takes 
place. Moreover we expect the correlation be- 
tween the two measurements of absorption to have 
a wide dispersion because the absorption mea- 
sured in the X-rays can be highly variable (see e.g. 
Risaliti et al.ll2007n . and because the two columns, 
in the X-ray and in the IR, are measured along dif- 
ferent paths as the emitting sources have different 
spatial locations and distributions. It is also quite 
plausible that, in addition to the absorption pro- 
duced by the medium that is responsible for the 
extinction seen in the IR, the X-ray source suffers 
also absorption by gas clouds, possibly ionized, 
close to the nucleus that may be dust-free, and 
thus not contribute to any emission or absorption 
in the IR. 

A s imilar comparison was presented bv lShi et al. 
( 20061 ) for a large compilation of AGNs of different 
type. Instead of the apparent Si optical depth, 
Tsi, they define the Si strength as an indicator 
of IR extinction. The Si strength is equivalent to 
e""^^'— 1, thus for our sources it ranges from —0.3 
to —0.8. We report their best linear fit obtained 
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for their entire AGN samples in Figure [14] after 
converting the Si strength to the apparent rsi- 
By comparing our sources with their relationship 
and their Fig. 3, we find that our sources oc- 
cupy the lower boundary of the Si strength-Nn 
relation they observe, and overlap with the most 
IR obscured sources in their sample which are all 
Seyfert 2 galaxies. Thus, the Nh and rsi of our 
sources are consistent with those found in other 
AGNs, but they overlap with those observed in 
the most IR obscured ones. This is a consequence 
of our sample selection. 

8.2. Dust covering fraction 

The dust covering fraction can be estimated by 
comparing the thermal reprocessed emission, e.g. 
L(6/im), and the luminosity of the heating source. 
The best proxy for the intensity of the heating 
source is given by the ultraviolet (UV) luminos- 
ity, since UV photons are more efficiently ther- 
malized than those at other wavelengths (both X- 
ray and IR). Since we do not know the intrinsic 
(before absorption) UV luminosity because it is 
mostly absorbed and only few sources are detected 
in the optical, we adopt the absorption-corrected 
X-ray luminosity as a proxy for the intensity of the 
heating source. Since the ratio between the UV 
and X-ray luminos ity varies with luminosity (e.g. 



Steffen et al. 20061 ) . in the sense that more lumi- 



nous AGN have lower X-ray/UV luminosity ratios, 
any trend of this ratio with the X-ray luminosity 
must be carefully interpreted. With this caveat 
in mind, we adopt the ratio between L(6/im) and 
Lo°3-8fcey ^ approximate indicator of dust cov- 
ering fraction. 

In Figure [151 we show the MIR over X-ray lumi- 
nosity ratios of our sample. In case of non X-ray 
detections, we show the X-ray luminosity corre- 
sponding to the upper limit to the flux. Since the 
flux is lower than the assumed upper limit, the in- 
trinsic luminosity could be much lower than esti- 
mated, but since we do not correct for absorption, 
it can also be higher than estimated. The lumi- 
nosity ratios for these sources are thus reported 
with an uncertainty of a factor of 10. 

We also report for comparison the luminosity 
ratios typical of AGNs (region enclosed by the 
dashed lines in Figure fTSl) and of the X-ray selected 
type 2 QSOs from Sturm ct al. (200f|) (hereinafter 
S06 sample). The range of ratios for AGNs was 
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Fig. 15. — Mid-IR (6/im rest-frame) over X-ray 
luminosity as a function of X-ray luminosity for 
the sources presented in this work (black full cir- 
cl es and open diamo nds) and the type 2 QSOs 



Sturm et al.l ( 2006f ) (black triangles). Full cir- 



cles are sources fitted with the T model and open 
diamonds are sources fitted with the T+C model. 
The X-ray luminosities have been corrected for ab- 
sorption. Large error bars on the luminosity ratio 
(a factor of 10 in the X-ray luminosity) are re- 
ported for sources for which only an upper limit 
to the X-ray flux is available. The dashed lines 
enclose the regio n where the rnajori ty of AGNs 
lie (adapted from jPolletta et al.ir2007t ). 
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obtained from a compilation of ratios for variou s 
types of AGN (adapted from lPolletta et al.ll2007h . 
These authors find that the observed range of 
MIR-to-X-ray luminosity ratios in AGNs is char- 
acterized by a wide dispersion, from 0.3 to 13, with 
no dependency on AGN type and luminosity, con- 
firming previous r e sults ( see e.g. iLutz et aLlf2004i : 



Sturm et al.ll2005L l2006t iHorst et al.ll2006h . The 



sources in the S06 sample are heavily absorbed in 
the X-rays and their MIR emission is mostly due 
to the AGN, but they do not show the Si feature 
in absorption. The ratios for the S06 sample are 
consistent with those typical of AGNs. 

Interestingly, most (9 out of 11) of the sources 
in our sample show higher ratios than the litera- 
ture samples and the S06 sample. The observed 
offset is not due to a different range of X-ray 
luminosities, as it is observed across the entire lu- 
minosity range of our sample. Yet, it is difhcult to 
establish whether there is a clear difference among 
our sample and the others. For example, the X- 
ray luminosity of LH_02 is likely underestimated. 
Indeed its soft X-ray spectrum and low X-ray lu- 
minosity suggest that its observed X-ray emission 
is scattered light and the intrinsic luminosity is 
completely absorbed. Moreover, the uncertainties 
associated with the absorption corrected X-ray 
luminosities can be quite large since these are 
typically faint sources and the column densities 
are derived from the hardness ratios rather than 
from X-ray spectral modeling. Deeper X-ray data 
for these sources or X-ray observations of similar 
depth as those already available for the rest of the 
sample would be necessary to confirm such a re- 
sult. If confirmed, this offset in Laz/jj/Lx might 
be simply due to a selection bias in our sample. A 
bias could have been introduced by the powerful 
AGN emission at MIR wavelengths which might 
be an indicator for a specific dust geometry, e.g. 
a more compact cloud distribution in the vicinity 
of the nucleus. Also the presence of Si absorption 
might have introduced a bias in favor of dustier 
sources and thus with larger dust covering fac- 
tors. However, we do not expect a correlation 
between deep Si features in absorption and cover- 
ing factors, i.e. large L(MIR)/L(X-ray) luminosity 
ratios, if the medium producing S i absorption i s 
clump y as sugge sted by Nenkov a et al.l (12002 ]: 



DuUemond fc van Bemmell (|2005h : iHonig et al 



not intercept any cloud, even if the covering frac- 
tion is large, and thus not exhibit the Si feature in 
absorption, whereas a source with a low covering 
factor might be seen through an optically thick 
cloud and show a deep Si feature. 

9. Discussion 

9.1. Starburst signatures in obscured AGNs 

For assigning the ultimate luminosity source of 
an obscured galaxy, it is necessary to know what 
fraction arises from an AGN and what fraction 
from a starburst. Our modeling discussed above 
indicates that the observed NIR and MIR emis- 
sions are mainly produced by dust heated by an 
AGN. Because of suggestions that the obscuring 
material may be associated with starbursts within 
the clouds that surround the AGN, we search for 
signatures of a starburst component in the spec- 
tra of these obscured AGNs. Because of poor S/N 
within the individual spectra (see Figure (2]), there 
are no confident detections of the features which 
are diagnostics of a starburst, such as PAHs and 
low ionization emission lines. Some of these fea- 
tures might be present in some spectra, but higher 
S/N spectra are required to confirm them. Here, 
we analyze the composite spectra in Figure [8] to 
search for weak PAH features that are not visible 
in individual objects. 

In none of the composite spectra, however, 
there are any confident detections of PAH fea- 
tures. Limits can be put on the strength of 
these features relative to the underlying contin- 
uum. The spectrum in Figure [8] with the high- 
est limit for a PAH feature is the composite for 
sources fit with the T model (panel a in Fig- 
ure [8]), in which a feature near 7.7/im might be 
a PAH feature whose peak F(7.7/im) exceeds by 
< 10% the flux density of the underlying con- 
tinuum at 7.7^m, F(7.7/im)=°"*. If this feature 
is real, the ratio F(7.7^m)/F(7.7/xm)™"*, where 
F(7.7/im) represents the flux of the total (contin- 
uum+line) at 7.7/im, would be 1.1. For objects 
classified as pure starbursts without any AGN 
contribution, t his ratio always exceeds 5 and is 



typica,llv ^10 (iGenze 



2006t IWeedman etall 



et al. 



2006a 



19981: iBrandl et al. 
Houck et al.ll2007l ). 



(|200fi[ ): IShi et al.l (|2nn(ih . Indeed, the LOS might 



Lower values are obse rved in some SMGs, but al- 
ways greater than 2 ijValiante et al.l 120071 ) . If a 



starburst is present in this spectrum, therefore. 
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it is very diluted by the underlying continuum 
of the AGN component, as seen in other AGNs 
where PAH emission appears to be weak or ab- 
sent in the p resence of a strong AGN continuum 
in the MIR dLutz et all Il998t iciavel et al] I2OO0I: 



Valiante et al 



2007t) . 



For use with spectra of poor S/N for which 
the total PAH fluxes are difficult to deter- 



mme. 



Houck et al.l ()2007f ) derive an empirical rela- 



tionship between the luminosity of the 7.7/^m fea- 
ture and the bolometric luminosity of a starburst, 
L f,S , using the s tarbu rst spectra and luminosities 



Brandl et all (|2006h . This relation is 



logL|^l^log[vLl''{^.^^lm)] 



0.78, 



(2) 



for vljf,^ (7.7/im) measured from the F(7.7/im)'™'^ 
at the peak of the 7.7/im feature. If the feature 
has lower equivalent width than in a pure star- 
burst, as in our composite spectrum, a correction 
is necessary to subtract the underlying continuum, 
F(7.7/L(m)'^°"*, from F(7.7/im) in order to deter- 
mine F(7.7/im)'™% and thus vhf,^ (7.7/im) that 
is due to the starburst. 

The composite spectrum in panel a of Figure [5] 
also indicates that j/F(7.7/im)^''"* — i^F(6/im). 
Therefore, the starburst 7.7/im luminosity can be 
derived as vLf,^ (7.7/im) < 0.1 xLcont{7-7fJ,m) 
= 0.1xL(6/im). Figure [3 and Table [3] illustrate 
that the median log{L{6iJ.m)/ergs^^) — 46.10, 
yielding a median log {i^L^^ (7.7 fim)) — 45.10 for 
the most luminous PAH feature that would sat- 
isfy the observed limit in the composite spec- 
trum. This limiting starburst would then have 
logL^^i/ergs'^ = 45.88 (equal to ^2xW^hQ or 
to a SFR~350 Mq yr~^) arising from the star- 
burst, using equation[2] This starburst luminosity, 
compared to the median AGN bolometric lumi- 
nosity, log L-^^^ /ergs^^ = 46.47, (see Figure [7] 
and Table [3|), indicates, therefore, that any star- 
burst which is present in the limit of the observed 
composite spectrum has, on average, a bolomet- 
ric luminosity < 26% of the bolometric luminosity 
arising from the AGN, and thus <20% of the total 
(AGN-|-Starburst) bolometric luminosity. 

Eight sources (out of 13 SWIRE sources and 
5 E-FLS sources) are detected in the FIR (70 
and/or 160/im). For these sources, the FIR lu- 
minosities (see Table [3|) are consistent with those 
observed in the most luminous starburst galax- 
ies ('--^lO^^Lo), but exceed the limiting starburst 



luminosities implied by the absence of PAH fea- 
tures (= lO*^-**ergss~-^=2xlO^^L0). In two cases 
(Nl_09, N2_06), these powerful FIR luminosities 
are consistent with those predicted by the torus 
models. However, in the other 6 cases there is 
evidence for an additional FIR component whose 
origin is still unclear. The FIR luminosity could 
arise from a starburst, or from dust far from the 
torus but still heated by the AGN. 

9.2. Obscuration at high-luminosities 

A comparison between the MIR properties of 
our sample with other samples from the literature 
shows that the selected sources can be as luminous 
as the most luminous MIR selected type 1 AGNs 
in the same redshift range. The estimated ratio 
between obscured and unobscured AGN at these 
luminosities (Lg^m ~10'*^ '^~'*^-^ ergs s~^) is about 
1.5-1.9, and the surface density is 17-22 obscured 
QSOs deg-2 to the limits of the SWIRE survey. 

However, as discussed in the previous sections, 
the ratio between obscured and unobscured QSOs 
derived here does not correspond to the ratio be- 
tween QSOs seen through the obscuring torus and 
those seen through the torus opening cone. In- 
deed, from our study, it emerges that the obscur- 
ing matter along the LOS is not always associated 
with the torus. Only about half (12 sources) of 
obscured QSOs are seen through the torus. Over- 
all and including the unobscured AGN popula- 
tion, we conclude that 35-41% of MIR selected 
red QSOs are unobscured, 37-40% are obscured 
by the torus, and the remaining 23-25% are ob- 
scured, but not by the torus. These values con- 
strain the covering factor of the torus at these ex- 
treme luminosities to be ~37-40%, and the half 
torus opening angle, measured from the torus axis, 
to be ~67°. This is significantly larger than the 
value of 46° estimated for F IR selected Seyfert 
galaxies ISchmitt et al.l (|200l[ ). Our results thus 
support the receding torus scenario, but are based 
on the interpretation that a large fraction of ob- 
scured QSOs are obscured by dust detached from 
the torus located along or outside the torus open- 
ing cone. 

The presence of foreground cold dust detached 
from the torus that might absorb the n uclear emis- 
sion was first suggested bvlKeell (|l98fl[ ) and, subse- 
quently supported by other stu dies at optical, IR, 



radio, and X-ray wavelengths ([Lawrence fc Elvid 
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1982; 



2006; 



Rigbv et al .120061: 1 Ogle et al.ll2 006':'Ma rtmez-Sansigrd;(Btisd roduce the observed NIR emission. The 



Brand et al.l2007HSaiina et al.h2007b : iUrrutia et al. 



20071 ). In our Galaxy, a gas-rich spiral, rsi is 
estim ated to be ^ 2. 5 in an equatorial direc- 
tion ( Rieke et al.lll989l ). It is possible, therefore, 
to produce the observed obscuration with dust 
in a normal spiral galaxy if this host galaxy is 
viewed edge-on, or in a very disturbed system as 
a merging galaxy. 

Whether the obscuring dust is in a gas-rich spi- 
ral host or in a starburst component, it clearly 
suggests that these extremely luminous and ob- 
scured AGNs are not hosted by dust-free elliptical 
galaxies as commonly s een in luminous AG Ns in 
the nearby universe fe.g. lDunlop et al1l2003[) . The 
sources in this study are thus powerful AGNs and 
they might be hosted by dusty galaxies. 

The cold absorber will also absorb the nuclear 
X-ray emission if located along its path. Thus, 
we can expect that all obscured QSOs, including 
those obscured only by the cold absorber and those 
obscured by the torus, will have large column den- 
sities in the X-ray (Nh>10^^ cm~^). Are our pre- 
dictions consistent with the findings from X-ray 
studies? 

The lack of X-ray data for the majority of 
the sources in our sample does not allow us to 
verify that these sources are also absorbed in 
the X-rays. Only 5 sources have available X- 
ray data, and they all show large column densi- 
ties (Nh>10^'^ cm~^), strongly suggesting that the 
same might be true for all these sources. At the 
X-ray luminosities of our sources, ~10^^ ergss"^, 
the observed ratio between X-ray absorbed and 
unabsorbed AGN in X-ray selected AGN sam- 
ples is, unfortun ately, still largely unconstrained. 



0.6-1.5:1 (e.g. lAkvlas et aiTl2006t iTaier et al 



20071 ). mainly because of the difficulty of detect- 
ing absorbed sources in the X-rays. These val- 
ues are lower than our predicted 1.5-1.9:1 ratio. 
On the other hand, when the selection effects 
against absorbed sources are minimized and the 
spectroscopic completeness is maximized, the ab- 
sorbe d to unabsorbed f raction can be as high as 
2.5:1 (jWang et al.ll2007l ). Since this ratio is higher 
than our estimates, it is more in agreement with 
our predictions. 

Instead of adding the cold absorber, an al- 
ternative way to reproduce the observed SEDs 
is to include an additional dusty component 



spectrum of this NIR component must decline 
steeply above ^Syitm, and thus have a very nar- 
row range of temperatures > 700 K. This NIR 
component could be produced by direct emis- 
sion from warm dust around the AGN and not 
be associate d with the torus, e.g. in the narrow 
line regi on tefstathi ou &: Rowan- Robinson 119951 : 
iRuiz et a l. 2000; Axon|l2Q0l|), but models do not 
usually reprodu ce the high tempera tures observed 
in our sample ( Groves et al. 20061 ). A more vi- 



able origin for this hot dust component might 
be scattered light in the NIR. This last hypoth- 
esis requires the presence of optically thin dust 
in the NLR which scatters the nuclear NIR emis- 
sion into the LOS (lEfstathiou fc Rowan-Robinson 



1995; Pier fc Krolik 



1993h . 



This scenario, ini- 



tially suggested theoreticallv (IPier fc Krolik 119931 : 
Efstathiou fc Rowan-RobinsoiJ 19951 ) . is also sup- 



ported by speckle interferometric NI R observa- 



tions of the Seyfert 2 galaxy NGC 1068 (jWeigelt et al 



20041 ). Reconstructed K-hand images show that 



the nuclear emission is extended parallel to the 
outflow and jet direction, instead of perpendicular 
as expected for the torus. Observations in polar- 
ized light would be necessary to test this scenario 
for our objects. Other possible explanations are a 
more complex geometry, or a different extinction 
curve. 

10. Summary and Conclusions 

In this work, we modeled the IR SEDs and 
IRS spectra of a sample of 21 MIR selected red 
AGNs at z=l. 3-3.0 with extreme MIR luminosi- 
ties {i'L^{6fim) = L{6^m) ~ 10**^ ergs s^^) with 
the goal of investigating the properties of their 
obscuring matter and dependence on luminosity. 
The sample was drawn from the three largest 
Spitzer surveys (SWIRE, NDWFS, fc FLS) by 
means of their extremely red infrared colors and 
presence of a Si absorption feature in their IRS 
spectra. Eighteen sources clearly show the Si ab- 
sorption feature with various optical depths (Fig- 
ure [2|). In the remaining three sources, the Si 
absorption is only tentatively observed. 

The obscured sources are as MIR luminous as 
the most luminous MIR selected type 1 AGNs cur- 
rently known at similar redshifts, but are charac- 
terized by Si optical depths typical of ULIRGs and 
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Seyfert 2 galaxies, that are usually less luminous 
MIR sources. The presented sample thus covers a 
new parameter space in the L(MIR)-TSi region. 

The observed SEDs and spectra are fit with 
the clu mpy torus models developed bv lHonig et al. 
( 20061 ). From the modeling, it emerges that the 



obscuring matter along the LOS is not always 
associated with the torus. Only about half (12 
sources) of obscured QSOs are seen through the 
torus and nine sources require, in addition to a 
torus, a cold absorber to reproduce the depth 
of the Si feature. The additional component is 
characterized by large optical depths (Ty'^=4-25). 
Seven of the sources fit with the additional cold 
absorber are characterized by tori at low inclina- 
tion angles, 9 =0-30°. The low inclination angles 
are required by the observed NIR emission which 
implies a view of the hottest dust components 
from the inner parts of the torus. Even though 
all sources show the Si feature in absorption, we 
conclude from the modelling that there is no pre- 
ferred torus inclination in our sample such that 
these obscured sources have been selected because 
they are observed through the torus. 

The best-fit models indicate that a compact 
non-flaring torus is preferred by the majority of 
the sources. The torus emission region is compact 
(more concentrated towards the nucleus) with the 
NIR emission being mainly due to dust in the 
vicinity of the nucleus. The preferred non-flaring 
torus in the majority of these luminous sources is 
also consistent wi th th e predictions of th e receding 
torus models (e.g. [Simpson 2005. : Honig fc Beckerti 
20071 ). 

Based on an estimate of the surface density 
of unobscured QSOs derived from MIR lumi- 
nosity functions ( Brown et al. 20061 ). and of ob- 
scured QSOs derived from a color-selected sam- 
ple of IR sources, the estimated surface den- 
sities are about 17-22 deg~^ for the obscured 
QSOs, and 11.7 deg~^ for the unobscured ones 
at z=l. 3-3.0. The estimated ratio between ob- 
scured and unobscured AGNs at these luminosities 
{L{6iim) ^lO'*^-^"^^-^ ergss"^) and redshift range 
is, thus, -1.5-1.9. OveraU we find that -35-41% 
of MIR selected red AGNs are unobscured, 23- 
25% are obscured but not by the torus, and the 
remaining 37-40% by the torus. The ratio of ob- 
scured to unobscured QSQs constrains the torus 
half opening angle, measured from the torus axis. 



to be —67°. This value is significantly larger than 
found for FIR select ed samples of A GN at lower 
luminosity (—46° Sc hmitt et al.ll200ll ). supporting 
the receding torus scenario. 

Five sources are detected in the X-rays, and 
large column densities (Nh>10^^ cm~^) are esti- 
mated for all of them. For these objects, the dust 
covering fraction estimated from the ratio between 
the MIR and X-ray luminosity are on average 
higher than what is usually measured in AGNs. 
However, the large uncertainties and limited sub- 
sample of detected X-ray sources, do not allow us 
to confirm such a result. Whether these sources 
are characterized by an excess of MIR emission 
with respect to the X-ray emission and whether 
this is due to a peculiar dust distribution with re- 
spect to the typical AGN population introduced 
by our selection is thus not clear. 

From the SEDs of the models, we estimate that 
the AGN bolometric luminosities, L^,^^ , range 
from 10^^ to 10^^ergss-\ and the AGN bolomet- 
ric luminosity can be derived from the 6/im lumi- 
nosity using the relation Log{L^^^)=Log{L{&^m))- 
0.32 ± 0.06. The bolometric luminosity of the 
whole system can be much higher, however, if the 
luminosity escaping the torus is taken into account 
and if an additional FIR component is included. 
Such a FIR component is not included in the AGN 
torus models, but it is observed in 8 cases. In half 
of the 8 FIR-detected sources the torus contribu- 
tion to the total IR luminosity is lower than that 
of the additional FIR component. For two of the 
remaining objects, the two contributions are simi- 
lar and in the other two, the torus contribution is 
larger. 

The estimated FIR luminosities for these 8 
sources are all greater than 3.3x lO^^ Lfri, which 
woul d imply SFRs of 600-3000 Mq yr"! (|Kennicutt 



19981 ) if these FIR luminosities arise from star- 
bursts. However, none of the sources show the 
PAH features characteristic of starbursts; a limit 
to the starburst component is estimated from the 
limit on PAH 7.7/im luminosity in the composite 
IR spectra. From this limit, we estimate that, 
on average, the contribution to the bolometric lu- 
minosity from a starburst characterized by PAH 
emission is <20% of the total bolometric lumi- 
nosity. This leaves unanswered the question of 
whether the extreme FIR luminosities of these 
sources arise from heavily obscured starbursts or 



31 



arise from AGN luminosity heating cool dust that 
is distant from the torus. 
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Table 1 

Optical and infrared data of the obscured QSOs sample 



CO 

00 



Survey 


Source 


iau 


Fu° 


(MJy) 








Fs.e^tm 


(MJy) 


Fs.S^m 

(MJy) 


Fs.Ofim 




Fyo^tin 


Fieo^im 


Rcf. 


Name 


Name 


Name 


(MJy) 


(MJy) 


(MJy) 


(MJy) 


(MJy) 


(MJy) 


(MJy) 


(mJy) 


(mJy) 




SWIRE 


LH.Ol 


SWIRE J104148. 93+592233.0 


<0.441 


<0.326 


<0.486 


<0.821 




5.4 


12 


<43 


120 


1351 


<17.0 


<90 


1 


SWIRE 


LH.02 


SWIRE J104605. 56+583742. 4 


<0.441 


<0.326 


<0.486 


<0.821 




12 


20 


69 


124 


1407 


<17.0 


<90 


1 


SWIRE 


LH_A4 


SWIRE J104409. 95+585224. 8 


0.360 


1.038 


1.186 


1.312 


1.823 


64 


150 


391 


1093 


4134 


<17.0 


<90 


2 


SWIRE 


LH_A5 


SWIRE J104453. 07+585453.1 


0.236 


0.333 


0.409 


0.611 


<1.000 


69 


141 


292 


467 


1190 


<17.0 


<90 


2 


SWIRE 


LH_A6 


SWIRE J104613. 48+585941. 4 


1.115 


1.799 


1.993 


2.850 


3.535 


27 


34 


58 


127 


1147 


<17.0 


<90 


2 


SWIRE 


LH.A8 


SWIRE J104528. 29+591326. 7 


<0.441 


0.517 


1.042 


2.223 


3.654 


32 


45 


91 


207 


2462 


9.8 


<90 


2 


SWIRE 


LH.A11 


SWIRE J104314. 93+585606. 3 


<0.411 


<2.473 


<2.324 


<3.269 




9.2 


22 


60 


117 


968 


<17.0 


<90 


2 


SWIRE 


N1.09 


SWIRE J160532. 69+535226. 4 


<0.801 


1.107 


1.491 


1.914 


<3.784 


15 


21 


<43 


73 


2894 


8.5 


<90 


1 


SWIRE 


N2.06 


SWIRE J163511. 43+412256. 8 


<0.801 


<0.409 


<0.752 


<1.324 


<3.784 


14 


18 


36 


116 


4231 


10.3 


<90 


1 


SWIRE 


N2.08 


SWIRE J164216. 93+410127. 8*= 


<0.801 


<0.409 


<0.752 


<1.324 


<3.784 


72 


176 


475 


1178 


3984 


<17.0 


80.5 


1 


SWIRE 


N2.09 


SWIRE J164401. 40+405715.0 


<0.801 


<0.409 


<0.752 


<1.324 


<3.784 


39 


58 


86 


143 


3715 


27.5 


121 


1 


NDWFS 


Bl 


SST24 ,1143001.91+334538.4" 




0.142 


0.302 


0.521 




8 


29 


102 


586 


3800 


<30.0 




1 


NDWFS 


B2 


SST24 ,1143644.22+350627.4'' 




0.325 


0.525 


0.753 




43 


97 


380 


732 


2300 


<30.0 




1 


NDWFS 


B3 


SST24 J142827. 19+354127. 7 




0.980 


2.290 


5.210 




296 


621 


1601 


3939 


10550 






1 


NDWFS 


B4 


SST24 J142648. 90+332927. 2° 




0.514 


0.912 


1.194 




53 


178 


493 


911 


2300 


<30.0 




1 


NDWFS 


B5 


SST24 J143508. 49+334739. 8 




0.390 


0.692 


0.905 




12 


21 


34 


237 


2600 


<30.0 




1 


FLS 


MIPS42 


SST24 J171758. 44+592816. 8 






0.186 






<9 


30 


103 


680 


4712 


12.8 




3 


FLS 


MIPS78 


SST24 J171538. 18+592540.1 






0.186 






<9 


39 


<72 


268 


2973 


<3.9 




3 


FLS 


MIPS15840 


SST24 J171922. 40+600500. 4 






0.462 






18 


25 


62 


197 


1821 


<4.8 




3 


FLS 


MIPS22204 


SST24 ,1171844.38+592000.5 






2.996 






<27 


39 


132 


478 


4101 


13.2 




3 


FLS 


MIPS22303 


SST24 J171848. 80+585115.1 






0.395 






<9 


<12 


<54 


103 


2030 


7.2 




3 



Note. — Typical uncertainties to the optical and IR fluxes are ~ 4%, and 10% of the measured fluxes, respectively. Upper limits correspond to 5ct. 
"The optical data correspond to Ug'r'i' z in SWIRE, to UBRI in NDWFS, and R in FLS. 
''Also known as ELAISC15 J164216. 9+410128. 
"N. 9 in lHouck et al.l 112003 1. 

''Also known as NDWFS J143644. 3+350627 and N. 13 in lHouck eFall I I2005I) . 
"Also known as FIRST ,1142648.9+332927 and ELAISR142650+332940B. 



References. — (1) this work; (2') IWeedman et al.l Il2006al) : (3) IRAC data from lSaiina et al.l I l2007al) . and R-band and MIPS[24] data from lYan et al.l II2007I) . 



Table 2 

irs low resolt'tion spectrosropy observation log 



Source 


Program 


SLl 


LL2 


LLl 


Name 


ID 


(seconds) 


(seconds) 


(seconds) 


LH.Ol 


136 


3x60 


6x120 


6x120 


LH.02 


136 


3x60 


6x120 


6x120 


LH.A4°' 


15 


2x60 


3x120 


3x120 


LH.AS" 


15 


3x60 


6x120 


6x120 


LH.Ae" 


15 


3x60 


6x120 


6x120 


LH.A8'' 


15 


2x60 


3x120 


3x120 


LH^ll* 


15 


5x60 


6x120 


6x120 


N1.09 


15 


3x60 


3x120 


3x120 


N2.06 


15 


2x60 


2x120 


2x120 


N2.08 


15 


2x60 


2x120 


2x120 


N2_09 


15 


2x60 


2x120 


2x120 


Bl 


12 


1x240 


2x120 


2x120 


B2 


15 


3x60 


3x120 


3x120 


B3 


15 


2x60 


2x120 


2x120 


B4 


15 


2x60 


3x120 


3x120 


B5 


15 


3x60 


3x120 


3x120 


MIPS42 


3748 


2x240 


2x120 


3x120 


MIPS78 


3748 


2x240 


2x120 


3x120 


MIPS15840 


3748 


2x240 


6x120 


7x120 


MIPS22204 


3748 


2x240 


2x120 


3x120 


MIPS22303 


3748 




2x120 


3x120 



Note. — IRS data processed with version 13.0 of the SSC pipeUne. 
^IRS data processed with version 11.0 of the SSC pipeline. 
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Table 3 

Infrared luminosities and optical depths 



o 



Source ID 


z 






is. 8 






r AGNa 


J AGN'b 


r SB C 


J AGN + 3 Bd 
Fl R 


r AGN + SBc 
^ [R 


J AGN+SBi 


fAGNs 
JlR 


fSBh 
JlR 




Tempi 


SFR"' 


LH.Ol 


2 


84 


46 


07 


46 


09 


46 


09 


46 


39 


46.45 


















>0.79 


>0.37 




LH_02 


2 


12 


45 


63 


45 


66 


45 


67 


45 


95 


46.00 


















1.98 


0.91 




LH_A4 


2 


54" 


46 


45 


46 


45 


46 


45 


46 


69 


46.83 


















0.87 


0.57 




LH_A5 


1 


94° 


45 


65 


45 


65 


45 


64 


45 


89 


46.08 


















0.63 


0.41 




LH_A6 


2 


20° 


45 


67 


45 


68 


45 


68 


45 


92 


46.02 


















0.43 


0.18 




LH_A8 


2 


42° 


46 


12 


46 


15 


46 


16 


46 


42 


46.49 


46.34 


46.36 


46.83 


46.86 





39 





61 


0.86 


0.33 


994 


LH_A11 


2 


30° 


45 


70 


45 


72 


45 


73 


46 


01 


46.08 


















0.85 


0.38 




N1_09 


2 


75 


46 


35 


46 


38 


46 


40 


46 


69 


46.72 


46.10 


46.19 


46.86 


46.88 





67 





33 


>0.86 


>0.42 


572 


N2_06 


2 


96 


46 


60 


46 


61 


46 


61 


46 


89 


46.92 


46.16 


46.20 


47.02 


47.04 





74 





26 


>0.67 


>0.20 


657 


N2.08 


2 


40 


46 


42 


46 


42 


46 


42 


46 


59 


46.69 


46.80 


46.81 


47.20 


47.23 





25 





75 


1.78 


1.14 


2868 


N2.09 


1 


98 


46 


05 


46 


08 


46 


09 


46 


41 


46.46 


46.75 


46.76 


47.12 


47.14 





20 





80 


1.01 


0.15 


2556 


Bl 


2 


46 


46 


41 


46 


41 


46 


41 


46 


56 


46.62 


















2.36 


1.49 




B2 


1 


77 


45 


91 


45 


91 


45 


90 


46 


09 


46.20 


















0.88 


0.71 




B3 


1.293° 


46 


29 


46 


30 


46 


30 


46 


53 


46.60 


















0.43 


0.21 




B4 


1.82 


45 


78 


45 


78 


45 


78 


45 


98 


46.11 


















1.72 


1.21 




B5 


2.00 


46 


01 


46 


03 


46 


04 


46 


26 


46.30 


















2.43 


1.39 




MIPS42 


1 


90P 


46 


20 


46 


21 


46 


21 


46 


44 


46.48 


46.13 


46.16 


46.72 


46.75 





52 





48 


0.42 


0.13 


613 


MIPS78 


2.43P 


46 


21 


46 


22 


46 


23 


46 


52 


46.58 


















1.04 


0.58 




MIPS15840 


2.23P 


45 


96 


45 


97 


45 


97 


46 


21 


46.28 


















0.37 


0.10 




MIPS22204 


2.08 


46 


17 


46 


19 


46 


20 


46 


48 


46.55 


46.18 


46.22 


46.77 


46.81 





52 





48 


1.04 


0.58 


688 


MIPS22303 


2.50P 


46 


11 


46 


13 


46 


14 


46 


28 


46.31 


46.34 


46.36 


46.78 


46.80 





31 





69 


1.10 


0.80 


994 



Note. — All luminosities arc logarithm in unit of crgss^^. Luminosities for a starburst component are estimated only for sources detected at 70 or 160/im. 

^ Logarithm of the AGN IR luminosity obtained by integrating the torus model between 3 and lOOO/^m. 

^ Logarithm of the AGN bolomctric luminosity obtained by integrating the torus model between 0.1 and 1000/^m. 

Logarithm of the starburst component FIR luminosity between 42.5 and 122.5/im. 

Logarithm of the AGN and starburst FIR luminosities between 42.5 and 122.5/J-m. 
^ Logarithm of the AGN and starburst IR luminosities between 3 and 1000//m. 
* Logarithm of the AGN and starburst bolometric luminosities. 
^Fraction of AGN luminosity to the total IR luminosity. 
^Fraction of starburst luminosity to the total IR luminosity. 

^Silicate optical depth derived from ^?^{-Pg'7^„i /-f^g 7^1^^,) with -fg^^^ estimated extrapolating a power-law model fit to the data at A <7/J-m in the rest-frame. 

^Silicate optical depth derived from IniFg^j^^ / Fq^j^^) with Fq^j^^ estimated from a type 1 QSO template normalized at the observed 24/j-m flux and redshiftcd at the redshift 
of the source. 

™Star formation rate derived from L^^^ using the IKennicuttI (|1998l) relationship. 



^^Optical spectroscopic icdshift for LH_A4 from fPoUetta et al.l l]2006() and for B3 from [Pesai et ai.i (Unni). 
°H^5: z = 1.89, LH.A6: 2.10, LH^8: 2.31, LH.A11: 2.25 in lWeedman et al.l 1 120063) . 

PMIPS42: 2=1.95±0.07, MIPS78: 2=2.65±0.1, MIPS15840: z=2.3±0.1, MIPS22303: z=2.34±0.14 in lYan et al.l 120071) . 



Table 4 
Torus model parameters 



Source ID 


Modor 




b" 






j^LOS! 






Tcmpti 




LH.Ol 


T 


3.0 


1 


.5 


60 


12 


6.0 








.29 


0.01 


0.016 


LH.02 


T 


3.0 


1 


.0 


90 


26 


26.0 





1 


.00 


0.32 


5.085 


LH.A4 


T 


3.0 


1 


.0 


60 


12 


9.0 





1 


.00 


0.32 


4.638 




T+C 


2.0 


1 


.0 


15 


26 


0.0 


6 


0, 


.62 


0.35 


0.575 


LH.A5 


T 


1.5 


1 


.5 


60 


18 


17.5 





0, 


.83 


0.27 


6.665 




T+C 


2.0 


1 


.0 


30 


26 


0.2 


7 


0, 


.62 


0.36 


0.703 


LH.A6 


T 


2.0 


1 


.0 


45 


26 


5.6 





0, 


.32 


0.07 


0.050 


LH.A8 


T 


3.0 


1 


.0 


60 


9 


6.5 





0, 


.81 


0.26 


0.078 


LH.A11 


T 


3.0 


1 


.0 


45 


12 


4.6 





0, 


.75 


0.25 


0.533 


N1.09 


T 


3.0 


1 


.0 


60 


12 


9.0 





1 


.04 


0.29 


0.021 


N2.06 


T 


1.1 


1 


.5 


30 


13 


1.0 





0, 


.34 


0.06 


0.004 


N2.08 


T 


2.0 


1 


.0 


90 


26 


26.0 





1 


.05 


0.40 


10.120 




T+C 


2.0 


1 


.0 


30 


26 


0.2 


21 


2, 


.33 


1.35 


2.515 


N2.09 


T 


3.0 


1 


.5 


45 


8 


5.2 





0, 


.49 


0.04 


2.307 


Bl 


T 


2.0 


1 


.0 


45 


26 


5.6 





0, 


.37 


0.11 


1.435 




T+C 


3.0 


1 


.0 


15 


26 


0.0 


25 


2, 


.94 


1.64 


0.042 


B2 


T 


2.0 


1 


.0 


45 


26 


5.6 





0, 


.37 


0.12 


0.391 




T+C 


2.0 


1 


.0 


30 


26 


0.2 


18 


1 


.85 


1.07 


0.540 


B3 


T 


2.0 


1 


.0 


45 


26 


5.6 





0, 


.37 


0.11 


0.530 


B4 


T 


2.0 


1 


.0 


30 


26 


0.2 





0, 


.01 


0.00 


10.000 




T+C 


2.0 


1 


.0 


30 


26 


0.2 


16 


1 


.60 


0.92 


4.898 


B5 


T 


2.0 


1 


.0 


45 


26 


5.6 





0, 


.56 


0.14 


0.006 




T+C 


1.5 


1 


.5 


90 


6 


6.0 


15 


2, 


.02 


0.93 


0.295 


MIPS42 


T 


2.0 


1 


.0 


45 


26 


5.6 





0, 


.43 


0.14 


0.390 


MIPS78 


T 


3.0 


1 


.0 


60 


12 


4.6 





1 


.00 


0.31 


2.153 




T+C 


3.0 


1 


.0 


30 


9 


0.3 


10 


0, 


.93 


0.40 


0.003 


MIPS15840 


T 


2.0 


1 


.0 


45 


26 


5.6 





0, 


.32 


0.07 


0.155 


MIPS22204 


T 


3.0 


1 


.0 


45 


12 


4.6 





0, 


.75 


0.25 


2.957 


MIPS22303 


T 


2.0 


1. 


.0 


90 


26 


26.0 





1 


.06 


0.40 


2.438 




T+C 


2.0 


1. 


.0 


90 


26 


26.0 


4 


1 


.54 


0.68 


0.775 



^T: torus model, T+C: torus+cold absorber model. The parameters for the T+C model are 
given only when the torus only (T) model gives a poor fit (x^ >0.1) to the observed SED and 
spectrum, or the T+C model significantly improves the • 

"^Index of the torus density radial distribution, 7ir(r) oc r^"", where r is the torus radius. 

"^Index of the clouds vertical distribution, H cc , where H is the torus scale height. 

'^Inclination of the torus axis with respect to the line of sight in degrees. 

^Average number of clouds along an equatorial LOS. 

^Average number of clouds al ong the LOS, which app roximately corresponds to the optical 
depth of the torus in the visual <|Natta &: PanagialllQS^ . 

^Cold absorber (CA) optical depth in the visible. 

^Absorber (torus for model T and torus+cold absorber for model T+C) optical depth asso- 
ciated with the Si feature measured by extrapolating a power-law mode 1 {tqI") or by fitting a 
QSO template {j'si^^) (see text). 

^Reduced obtained from the best-fit model and the observed data at A >l/im in the 
rest-frame. 
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Table 5 
x-ray properties 



Source ID 






Log(Lo.3_sfeev) 






(10"^'^ orgscm^^g-i) 


(10^2 cm-2) 


(crgss^^) 


( orgs s ^ ) 


LH.02 


4.51 


12+"-^° 


44.07 


44.08 


LH.A4 


1.95 


105+^f 


43.83 


45.08 


LH.A5 


19.31 


9 4+*3 


44.61 


45.06 


LH.A6 


24.94 


6.3;|;= 


44.78 


45.10 


LH.A8 


4.45 


40+29 


44.19 


45.00 


LH.A11 


<1.0 


<43.46 




Bl 


<10.0 




<44.53 




B2 


<10.0 




<44.21 




B3 


<10.0 




<43.89 




B4 


<10.0 




<44.23 




B5 


<10.0 




<44.33 





^Broad-band (0.3-8 kcV) X-ray flux. 

^Eff'ectivc column density derived from the hardness ratios assuming that the intrinsic spectrum 
is a power-law model with photon index, F equal to 1.7 and Galactic Nh (6x10^^ cm~^). Uncer- 
tainties reflect only the statistical errors from the observed counts and do not include uncertainties 
in the spectral model. 

■^Absorption-corrected X-ray luminosities derived assuming the effective column density, Nh^'''"^ , 
and an intrinsic spectrum modeled with a power-law model with photon index, F— 1.7. 

Note. — For more details on the X-ray data see lPolletta et al.l l]2006|) . Upper limits to the X-ray 
fluxes correspond to 90% completeness. 
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